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ABSTRACT 
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In  an  earlier  AFOSR  funded  investigation,  experimental  research  was  conducted  to 
determine  the  influence  of  free-stream  turbulence  on  turbulent  boundary  layer  heat 
transfer  and  mean  profile  development.  The  data  obtained  under  this  earlier  contract 
indicated  that  both  the  skin  friction  and  the  heat  transfer  increased  significantly 
with  increased  free-stream  turbulence  level.  Under  the  present  investigation,  de¬ 
tailed  boundary  layer  turbulence  structural  data  and  turbulent  heat  transfer 
data  were  obtained  for  experimental  test  conditions  and  profile  locations  selected 
from  the  earlier  test  matrix.  Numerous  measurements  assured  that  the  present  test 
conditions  (boundary  layer  development  and  free-stream  turbulence  distributions) 
duplicated  those  of  the  earlieKAFOSR'contract.  The  purposes  for  making  these  present 
detailed  boundary  layer  turbulence  measurements  were:  (1)  to  provide  data  to  which 
current  finite-difference  boundary  layer  turbulence  models  could  be  compared,  and  (2) 
to  generate  a  data  base  for  the  development  of  new  analytical  models  for  boundary 
layer  heat  transfer  prediction.  .  The  results  from  the  present  program  have  shown  that 
the  distributions  of  both  the  turbulence  kinetic  energy  and  the  turbulence  structural 
coefficients  were  affected  by  increased  levels  of  free-stream  turbulence.  Local  pro¬ 
file  measurements  indicated  that  the  effect  of  increased  free-stream  turbulence  was 
to  decrease  the  near-wall  turbulent  Prandtl  number  relative  to  values  expected  for 
low  free-stream  turbulence.  Turbulent  Prandtl  numbers  in  the  outer  region  of  the 
boundary  layer  were  slightly  increased  for  higher  free-stream  turbulence.  A  turbu¬ 
lence  dependent  correlation  for  the  measured  distribution  of  turbulent  Prandtl  number 
is  given. 

With  the  completion  of  the  experimental  portion  of  this  investigation,  a  theo¬ 
retical  effort  was  made  to  assess  the  capability  of  a  finite  difference  boundary 
layer  computer  program,  ABLE  (Analysis  of  the  Boundary  Layer  Equations)  for  predict¬ 
ing  the  effect  of  free-stream  turbulence  on  momentum  and  thermal  boundary  layers.- 
Comparisons  with  experimental  data  of  mean  flow  velocity,  mean  flow  temperature, 
Reynolds  shear  stress,  turbulent  heat  transport,  and  turbulence  kinetic  energy  were 
made  in  this  investigation.  In  addition,  the  turbulent  Prandtl  number  correlation 
deduced  from  the  experimental  measurements  was  used  in  the  boundary  layer  analysis 
and  its  effect  on  surface  heating  evaluated.  The  results  indicated  that  this  boun¬ 
dary  layer  analysis,  which  uses  a  one  equation  eddy  viscosity  turbulence  model,  can 
provide  adequate  predictions  of  zero  pressure  gradient  flows  with  high  free-stream 
turbulence  and  wall  heating. 


INTRODUCTION 


The  search  for  improved  gas  turbine  performance  has  led  steadily  in  the  direc¬ 
tion  of  higher  turbine  inlet  temperatures.  The  last  twenty  years  have  seen  an  in¬ 
crease  in  turbine  inlet  temperatures  of  roughly  1400°F  but  an  increase  in  allowable 
blade  metal  temperature  of  only  roughly  200°F.  The  difference  between  these  two 
increases  in  temperature  can  be  related  directly  to  improved  cooling  technology.  As 
an  integral  part  of  this  advancing  cooling  technology,  engine  manufacturers  are  con¬ 
tinually  seeking  improved  techniques  for  calculating  heat  transfer  coefficient  dis¬ 
tributions  on  gas  turbine  airfoils.  As  the  level  of  cooling  technology  has  been 
driven  upward,  and  with  it  turbine  inlet  temperature,  it  is  not  surprising  that  the 
result  is  a  design  methodology  which  is  extremely  unforgiving  of  even  small  errors. 
The  temptation  is  always  present  to  overcool  the  airfoils  but  this  is  unacceptable 
due  to  the  powerful  negative  impact  of  cooling  air  on  the  cycle  and  on  turbine  effi¬ 
ciency.  It  is  this  dilemma  which  has  often  led  to  extremely  long  and  expensive  de¬ 
velopmental  testing  of  advanced  technology  turbines. 

Gas  turbine  thermal  design  systems  are  typically  not  based  on  fundamental  fluid 
mechanics  and  heat  transfer  data  and  analysis  alone  but  rather  they  are  calibrated, 
or  adjusted,  to  provide  agreement  with  engine  experience.  Without  the  benefit  of  a 
first-principles  understanding  of  the  effects  involved  there  is  the  likelihood  that 
a  designer  will  unknowingly  either  overcool  the  component  or  go  beyond  the  range  of 
validity  of  the  design  system  calibration.  There  is,  then,  a  clear  requirement  for 
the  development  of  airfoil  heat  transfer  distribution  prediction  procedures  which 
are  based  on  fundamental  fluid  mechanics  and  heat  transfer  data.  The  great  emphasis 
placed  on  the  development  of  accurate  boundary  layer  calculation  techniques  over  the 
past  few  years  reflects  the  recognition  of  these  needs. 

One  particularly  important  topic  in  the  general  context  of  turbine  airfoil  con¬ 
vective  heat  transfer  is  the  influence  of  the  free-stream  turbulence  on  fully  turbu¬ 
lent  boundary  layer  development.  It  has,  of  course,  long  been  recognized  that  in¬ 
creasing  the  free-stream  turbulence  level  can  cause  a  forward  shift  of  the  laminar 
to  turbulent  transition  region.  This  particular  phenomenon,  the  reduction  of  the 
boundary  layer  transition  Reynolds  number  with  increased  free-stream  turbulence 
level,  is  well  documented  in  the  open  literature  for  zero  pressure  gradient  flows 
and  can  be  adequately  predicted  with  currently  available  boundary  layer  prediction 
schemes.  In  addition,  a  number  of  investigators  have  studied  the  effects  of  free- 
stream  turbulence  level  on  turbulent  boundary  layer  growth,  profile  structure,  skin 
friction  distribution  and  heat  transfer.  The  consensus  of  these  studies,  is  that 
free-stream  turbulence  has  a  very  large  and  important  influence  on  both  the  heat 
transfer  and  the  boundary  layer  characteristics.  As  an  example,  it  has  been  shown 
in  a  recently  completed  AFOSR  funded  contract  at  UTRC  that  a  free-stream  turbulence 
intensity  of  5  percent  produces  an  increase  in  Stanton  number  of  approximately  15 
percent  over  the  value  expected  for  a  low  turbulence  freestream.  While  a  number  of 
existing  boundary  layer  analysis  procedures  (including  the  UTRC  ABLE  code)  account 
reasonably  well  for  the  influence  of  free-stream  turbulence  on  mean  velocity  profile 
development  and  skin  friction,  no  currently  available  analysis  satisfactorily  pre¬ 
dicts  the  observed  increased  heat  transfer  rates. 
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The  present  program  was  designed  to  provide  detailed  boundary  layer  ti  lulence 
and  turbulent  heat  flux  distribution  data  for  a  range  of  free-stream  turbulence 
levels.  As  part  of  this  program  these  experimental  data  were  employed  to  evaluate 
the  analytical  turbulence  models  currently  incorporated  in  the  UTRC  ABLE  code.  It 
is  anticipated  that  in  the  future  these  experimental  data  will  be  used  by  both  UTRC 
and  other  workers  in  the  field  of  boundary  layer  computation  for  development  of  new 
analytical  turbulence  models. 

The  contract  effort  consisted  of  acquiring,  documenting  and  analyzing  experi¬ 
mental  flat  wall  boundary  layer  mean  and  fluctuating  profile  data  to  determine  the 
influence  of  free-stream  turbulence  on  fully  turbulent  boundary  layer  flows.  For 
fully  turbulent,  zero  pressure  gradient  flows,  the  following  profile  data  were  ob¬ 
tained  for  a  range  of  free-stream  turbulence  intensities;  boundary  layer  mean  and 
fluctuating  velocities  and  temperatures,  turbulent  shear  stresses,  and  turbulent 
Prandtl  numbers.  In  addition,  in  order  to  improve  the  ability  of  the  UTRC  boundary 
layer  deck  to  predict  the  effects  of  free-stream  turbulence  on  heat  transfer  rates, 
a  turbulent  Prandtl  number  distribution  model  was  incorporated  into  the  UTRC  boun¬ 
dary  layer  code.  Calculations  were  carried  out  employing  the  measured  turbulent 
Prandtl  number  distributions  and  comparisons  made  between  the  predicted  and  measured 
heat  transfer  distributions. 
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DESCRIPTION  OF  TEST  EQUIPMENT 


I.  Wind  Tunnel  and  Heat  Transfer  Test  Surface 

All  experimental  data  for  the  present  investigation  were  obtained  in  the  United 
Technologies  Research  Center  (UTRC)  Boundary  Layer  Wind  Tunnel  (Fig.  1).  This  tunnel 
was  designed  specifically  to  generate  large-scale,  two-dimensional,  incompressible 
boundary  layers  with  Reynolds  numbers  and  free-stream  turbulence  levels  typical  of 
turbomachinery  airfoils.  Complete  descriptions  of  this  facility  including  measure¬ 
ments  documenting  the  tunnel  flow  uniformity  and  two-dimensionality  of  the  test 
boundary  layers  are  given  in  Refs.  1  and  2.  The  tunnel  test  section  consisted  of  a 
flat  upper  wall  instrumented  for  heat  transfer  measurements  which  served  as  the 
boundary  layer  test  surface,  plexiglass  vertical  sidewalls  and  a  flexible  lower  wall. 
The  test  section  was  34-in.  wide,  96-in.  long  and  8-in.  high  at  the  entrance.  For 
all  test  flows  in  this  study  the  lower  flexible  test  section  wall  was  adjusted  to 
produce  a  constant  velocity  along  the  test  section. 

A  photograph  of  the  Boundary  Layer  Wind  Tunnel  is  presented  in  Fig.  2.  Also 
shown  in  Fig.  2  are  both  the  telescope  used  to  position  probes  relative  to  the  test 
wall  and  the  computer  controlled  probe  traverse  mechanism. 

The  boundary  layer  test  surface  (upper  wall  of  test  section)  consisted  of  a  uni¬ 
form  heat  flux  electrically  heated  plate  instrumented  for  the  measurement  of  local 
convective  coefficients  (Fig.  3).  The  heated  flat  plate  was  constructed  from  a  block 
of  rigid  urethane  foam  34- in.  wide  by  96-in.  long  by  4-in.  thick  mounted  in  a  plexi¬ 
glass  frame  with  6-in.  wide  strips  of  316  stainless  steel  foil  cemented  to  the  test 
surface.  Details  of  the  flat  plate  model  and  its  instrumentation  are  presented  in 
Refs.  1  and  2.  Rigid  foam  was  employed  for  the  substrate  of  the  heated  flat  plate 
model  because  of  its  extremely  low  thermal  conductivity  (0.025  Btu/hr  ft°F).  Less 
than  1/2  percent  of  the  heat  generated  on  the  surface  of  the  plate  was  conducted  aw 
from  tne  test  surface.  Electric  current  passing  through  the  metal  foil  strips 
cemented  to  the  test  surface  produced  the  surface  heating.  The  metal  foil  strips 
were  wired  in  series  and  were  powered  by  a  single  low-ripple,  regulated  dc  power 
supply.  The  foil  test  surface  was  instrumented  with  an  array  of  203  Cr-Al  0.13 
mm  diameter  bead  welded  thermocouples.  Each  thermocouple  was  welded  to  the  back 
surface  of  the  foil  through  a  hole  in  the  rigid  foam  plate.  Forty-eight  surface 
static  pressure  taps  were  also  installed  along  the  test  surface. 

The  dc  power  current  passing  through  the  surface  strips  was  measured  using  two 
precision  shunt  resistors  and  a  digital  voltmeter.  The  temperatures  of  the  test  sur¬ 
face  thermocouples  were  measured  relative  to  a  single  test  section  free-stream  refe¬ 
rence  junction  using  a  digital  voltmeter. 

In  order  to  insure  a  known,  constant  test  surface  emissivity  and  hence  a  known 
radiation  loss,  the  completed  foil  test  surface  was  coated  with  3M  C-101  high  emissi¬ 
vity  flat  black  paint  (e  =  0.99).  Test  results  indicated  that  this  surface  was  aero- 
dynamically  smooth,  producing  no  premature  boundary  layer  transition. 
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Local  convective  coefficients  were  determined  by  ignoring  the  negligible  cr- 
duction  losses,  subtracting  power  lost  through  thermal  radiation  and  dividing  by  the 
temperature  difference  from  the  wall  (Tw)  to  the  freestream  (Te).  To  illustrate  the 
magnitude  of  the  radiation  losses  from  the  test  surface,  for  Ue  =  100  ft/sec,  for 
turbulent  boundary  layer  flow  with  Tw-Te  =  25°F,  the  radiation  loss  was  approximately 
4  percent  of  the  total  surface  power. 

As  shown  in  Fig.  1,  at  the  test  section  entrance  a  bleed  scoop  formed  the  lead¬ 
ing  edge  of  the  heated  boundary  layer  test  surface.  The  scoop,  which  was  mounted 
smoothly  on  the  front  edge  of  the  heated  test  wall,  provided  a  very  short  unheated 
starting  length  (c  =  1.7  in.)  upstream  of  the  heated  test  surface.  The  leading  edge 
of  the  scoop  was  a  4  x  1  ellipse  in  order  to  prevent  a  local  separation  bubble  and 
premature  transition  of  the  test  surface  boundary  layer.  Details  of  the  scoop  in¬ 
cluding  its  instrumentation  and  adjustment  are  given  in  Refs.  1  and  2. 


2.  Turbulence  Generating  Grids 

As  described  in  Ref.  1,  this  wind  tunnel  has  a  relatively  low  residual  test  sec¬ 
tion  turbulence  level  (<  1/4%).  Higher  turbulence  levels  required  for  this  study 
were  generated  by  inserting  various  square  array  biplane  grids  constructed  from  rec¬ 
tangular  bars  at  the  entrance  to  the  main  tunnel  contraction  (see  Fig.  1).  Four  tur¬ 
bulence  generating  grids  were  designed  using  the  correlations  of  Ref.  3.  The  grids 
will  be  referred  to  as  Grids  1,  2,  3  and  4  corresponding  to  mesh  widths,  M,  of  7/8, 

2  9/16,  7  and  9  in.  The  minimum  turbulence  configuration  (no  grid)  will  be  referred 
to  as  Grid  0.  Details  of  the  grid  configuration  are  given  in  Ref.  1.  This  present 
arrangement  differs  from  that  used  for  nearly  all  the  earlier  investigations  of  this 
subject  in  which  the  turbulence  grids  were  located  in  the  test  section  just  upstream 
of  the  boundary  layer  test  surface.  The  benefits  derived  from  locating  the  grids  at 
the  contraction  entrance  were  that  the  generated  turbulence  was  more  homogeneous  and 
had  a  lower  decay  rate  along  the  test  section.  Since  grid  generated  turbulence  decays 
approximately  as  u'/U  a  (x/b)"^77  (Ref.  3),  the  change  in  turbulence  level  with  dis¬ 
tance  along  the  test  section  was  reduced  by  increasing  the  distance  from  the  grid  to 
the  test  section  entrance.  In  addition,  the  results  of  Ref.  3  indicate  that  approxi¬ 
mately  10  grid  mesh  lengths  are  required  to  establish  a  uniform  turbulent  flow.  Lo¬ 
cating  the  grid  a  distance  upstream  of  the  test  section  requires,  of  course,  a  more 
coarse  grid  to  achieve  a  given  test  section  turbulence  intensity. 

Another  effect  considered  was  the  expected  influence  of  the  contraction  of  the 
components  of  the  grid  generated  turbulence.  It  was  recognized  that  rearrangement 
of  the  relative  magnitudes  of  the  turbulence  components  would  occur  due  to  the  con¬ 
traction.  However,  since  the  contraction  ratio  was  small  (2.8),  it  was  concluded 
that  any  effects  of  induced  anisotropy  would  be  small  in  comparison  to  the  advantages 
gained  in  homogeneity  and  reduced  decay  rate.  To  determine  the  validity  of  the 
assumption,  all  three  components  of  the  test  section  turbulence  were  documented  for 
all  test  cases. 


3.  Boundary  Layer  Total  Pressure  and 
Thermocouple  Probes  and  Traverse  Control 


Boundary  layer  mean  velocity  profile  data  were  measured  using  United  Sensor 
Model  Ba-0.020  impact  probes  with  flattened  tips.  The  probes  used  in  the  program 
were  inspected  for  defects  using  both  a  Nikon  Model  II  toolmakers  microscope  and  a 
Jones  and  Lamson  Model  PC14  Shadowgraph.  Mean  temperature  data  were  measured  with 
miniature  thermocouple  probes  designed  using  the  results  of  Ref.  4.  The  thermo¬ 
couple  sensing  elements  for  these  probes  were  constructed  from  0.001  in.  dia.  Chromel- 
Alumel  bead  welded  wires.  The  thermocouple  bead  ( *  0.003  in.  dia.)  was  located  at 
the  center  of  the  probe  support  prongs  which  were  fabricated  of  heavier  Chromel 
and  Alumel  wire.  The  results  of  Ref.  4  indicate  that  a  probe  of  this  design  was 
virtually  free  of  wire  conduction  errors  and  was  capable  of  measuring  boundary 
layer  mean  temperature  profile  data  into  the  viscous  sublayer  region. 

Movement  of  the  boundary  layer  probes  was  achieved  using  and  L.C.  Smith  ball/ 
screw  traverse  drive  with  an  optical  shaft  encoder  capable  of  resolving  relative  probe 
location  to  within  0.0005  in.  The  traverse  mechanism  was  suspended  on  a  linear  ball 
bearing  track  beneath  the  test  section.  The  traverse  could  be  located  anywhere  in 
the  center  75  percent  of  the  test  section  width  from  the  leading  to  trailing  edges  of 
the  test  wall.  A  telescope  sighted  through  the  tunnel  sidewall  was  used  to  accurately 
position  probes  relative  to  the  test  wall.  Estimated  absolute  accuracy  of  measured 
probe  distance  from  the  test  surface  was  0.002  in.  for  any  location  in  the  test  boun¬ 
dary  layers. 
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HOT  WIRE  DATA  ACQUISITION  AND  ANALYSIS  TECHNIQUES 


1.  General 

Measurements  of  fluctuating  velocities  and  temperatures  in  the  test  boundary 
layers  were  obtained  using  multi-element  hot  wire  anemometry  techniques.  For  a  large 
number  of  the  test  cases  the  wind  tunnel  was  operated  without  wall  heating,  the 
resulting  boundary  layers  being  isothermal.  Both  vertical  and  horizontal  x-type  2 
wire  probes  were  employed  for  these  isothermal  test  cases.  For  cases  with  wall 
heating  the  velocity  and  temperature  fluctuations  in  the  flows  were  determined  by 
using  specially  designed  3  wire  probes,  one  wire  of  which  was  operated  at  a  lower 
overheat  than  the  other  two.  Detailed  descriptions  of  both  the  2  and  3  wire  probes 
are  given  below  in  section  2.  The  voltage  signals  from  the  various  hot  wire  probes 
were  digitized,  recorded  and  subsequently  reduced  to  fluctuating  velocity  and  tem¬ 
perature  records  using  a  minicomputer.  A  detailed  description  of  the  data  system 
is  provided  in  section  4.  An  analysis  of  the  uncertainties  of  the  various  hot  wire 
measurements  is  given  in  Appendix  A. 


2.  Description  of  the  Hot  Wire  Probes 


2. 1  Probe  Design 

The  present  study  involved  the  measurement,  using  arrays  of  inclined  hot  wires, 
of  fluctuating  velocities  and  temperatures  within  boundary  layer  flows.  In  order  to 
minimize  potential  errors  for  these  measurements  (errors  largely  arising  from  the 
inherent  mean  velocity  and  temperature  gradients  in  the  flows  and  the  finite  probe 
size)  the  hot  wire  probes  were  custom-designed  and  fabricated  specifically  for  this 
program.  The  results  from  a  large  number  of  previous  boundary  layer  turbulence  and 
general  hot  wire  studies  were  incorporated  into  the  probe  designs  (Refs.  5-14).  For 
the  2  wire  x-type  probes  used  in  the  isothermal  tests  three  important  design  principles 
were  adopted  from  these  earlier  studies.  (1)  To  reduce  the  effects  of  the  mean 
gradients  in  the  flows  the  active  length  (or  the  size  of  the  array  of  wires  in  the 
direction  of  the  gradients)  of  the  wires  should  be  minimized.  (2)  To  reduce  end 
effects  (nonuniform  temperature  along  the  active  length)  and  to  insure  that  a  "Cham¬ 
pagne  k2"  (Ref.  7)  form  of  angular  sensitivity  could  be  employed,  an  active  length/ 
diameter  ratio  of  200  was  chosen.  (3)  To  maximize  the  spatial  correlation  coefficient 
(maximum  accuracy  of  cross-products  such  as  Reynolds  stress)  without  introducing  wire 
cross-talk  effects  a  transverse  wire  spacing  of  3z/4  was  chosen. 

Considerations  (1)  and  (2)  required  that  the  diameter  (d)  be  as  small  as  pos- 
$ible--the  limitation  being  practical  considerations  of  probe  fabrication  and  sensor 
survivability.  A  probe  development  program  (UTC  funded)  conducted  jointly  with  DISA, 
Inc.  resulted  in  the  conclusion  that  the  minimum  practical  wire  diameter  for  these 
probes  was  2.5  urn  (0.0001  in.)  for  platinum  plated  tungsten  wires.  From  considera¬ 
tion  (2)  the  active  length  of  the  wires  was  chosen  as  .50  mm  (0.020  in.)  and  from 
(3)  the  transverse  spacing  was  selected  as  0.015  in.  These  wire  arrays  were  employed 
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for  x-type  configurations  oriented  in  both  the  vertical  and  horizontal  planes.  As 
will  be  discussed  in  the  results  section,  cross-checks  indicate  that  the  fluctuating 
data  measured  with  probes  of  this  design  are  consistent  and  accurate. 

The  special  3  wire  probes  consisted  of  vertical  x-type  wire  arrays  with  a  third 
wire  mounted  equidistant  between  the  wires  of  the  x.  This  third  wire  was  parallel 
to  one  of  the  wires  of  the  x  array.  All  three  wires  were  constructed  from  the  same 
material  (platinum  plated  tunsten)  and  had  the  same  diameter  (0.0001  in.)  and  active 
length  (0.020  in.).  The  transverse  separation  between  adjacent  wires  of  the  3  wire 
array  was  0.015  in.  With  this  wire  arrangement  the  two  parallel  wires  of  the  array 
were  exposed  to  equal  effective  velocities  during  any  given  data  sample  period.  De¬ 
tails  of  the  techniques  used  to  determine  instantaneous  velocities  and  temperatures 
with  the  3  wire  probes  are  given  in  section  4. 

2.2  Probe  Calibration 

Prior  to  calibration  all  probe  sensors  were  operated  for  approximately  two  hours 
in  the  100  ft/sec  mainstream  of  the  wind  tunnel.  During  this  "wire  curing"  step  the 
sensors  were  set  to  operate  at  overheat  ratios  slightly  higher  than  those  used  during 
actual  testing.  These  "curing"  steps  (1)  provided  some  assurance  that  the  sensors  on 
a  given  probe  would  be  likely  to  survive  the  calibration  and  testing  environments, 
and  (2)  improved  the  stability  of  the  calibration  constants  of  a  given  sensor.  Each 
probe  was  calibrated  for  temperature-resistance  characteristics  in  a  low  temperature 
recirculating  oven.  Typically  five  temperature  vs.  resistance  points  were  measured 
for  each  sensor.  A  least-squares  data  reduction  program  was  used  to  find  a  best 
temperature-res i stance  coefficient . 

Rw:R32['  +  0(Tw-32)]  (1) 

where  Rw  =  resistance  of  the  active  sensor 
R32  *  sensor  resistance  at  32°F 
Tw  =  sensor  temperature 
a  =  temperature-resistance  coefficient 

Following  the  preliminary  "burn-in"  and  the  resistance  temperature  calibration 
each  sensor  was  calibrated  for  velocity  and  angular  sensitivity  in  a  low- turbulence 
H-in.  dia.  jet  flow.  The  sensors  of  the  2  wire  probes  were  calibrated  to  an  over¬ 
heat  (Rw-hot/Radiabatic)  1*5.  For  the  3  sensor  probes  the  outside  two  sensors 
were  calibrated  at  an  overheat  of  1.5  while  the  center  sensor  was  calibrated  at  an 
overheat  of  1.2.  With  the  main  probe  support  stem  oriented  perpendicular  to  the  jet 
axis  (wires  ±  45°  to  the  jet  axis)  mean  velocity  and  bridge  output  voltage  were  re¬ 
corded  for  approximately  20  jet  speeds  ranging  from  7  to  130  ft/sec.  The  mean 
response  equation  of  each  sensor  was  assumed  to  be  of  the  form 

Nu= A|+B|Re0,49  (2) 
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which  can  be  algebraically  manipulated  to 


A2tR,+Rw)  Q76 
R» 


tTw-T)  + 


B2(Rt+Rw)< 


(Tw-T)  U 


0.45 

E 


(3) 


where  Ew  =  wire  voltage 

Rs  *  probe  body,  cable  and  internal  anemometer  resistance 
Rw  *  sensor  resistance 
T  =  air  temperature 
Tw  =  sensor  temperature 
U£  =  effective  velocity 
A2,  B2  =  empirical  constants 

The  constants  A2  and  B2  were  determined  for  each  sensor  from  a  least-squares  fit 
of  the  data  to  Equation  (3).  Next,  using  a  pitching  fixture,  pitch  angle  versus  vol¬ 
tage  data  were  obtained  with  the  probes  rotated  from  +20°  to  -20°  in  steps  of  5°. 

The  center  of  pitch  coincided  with  the  intersection  of  the  wires  of  the  x.  Pitch 
sensitivity  data  were  obtained  for  three  jet  velocities,  50,  80  and  100  ft/sec.  The 
angular  sensitivity  of  the  wires  was  assumed  to  conform  to  Champagne’s  k?  law  (Ref.  7), 

UE2  ($)  =  U2  [<f>  0 )  (cos2  <f>  +  k2  sin2  $)  ( 4 ) 

where  $  =  angle  between  wire  and  direction  normal  to  the  flow  (;  ~  45° 
with  wall  probe  stem  normal  to  the  flow) 

*  effective  velocity 

Using  a  least-squares  routine  to  find  a  best  fit  of  the  pitch-voltage  data  to 
Champagne's  equation,  optimum  values  of  k  were  determined  for  each  sensor. 

In  summary,  the  temperature-resistance,  mean  velocity  and  pitch  calibrations 
were  used  to  determine  the  following  calibration  constants. 

(1)  R32  -  sensor  resistance  at  32°F 

(2)  a  -  temperature-resistance  coefficient 

(3)  A2  and  B2  -  empirical  constants  (Eq.  3) 

(4)  k  -  empirical  constant  (Eq.  4). 


3.  Description  of  the  Data  System 

For  all  test  cases,  both  for  isothermal  flows  and  for  flows  with  wall  heating, 
the  multi-element  hot-wires  were  driven  by  Thermo  Systems,  Inc.  (TSI)  Model  1050 
constant  temperature  anemometers.  Signals  from  the  anemometers  were  first  passed 
through  a  wide  band  amplifier  (Preston  Model  8300  XWB)  and  then  digitized  using  a 
TSI  Model  1075  Multichannel  Digitizer.  A  feature  of  this  particular  analog-to- 
digital  converter  which  is  important  to  this  application  is  that  the  various 
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channels  are  sampled  and  held  simultaneously.  This  simultaneous  sample-hold  feature 
permits  cross-products  of  the  various  fluctuating  quantities  to  be  computed.  Storage 
restrictions  of  the  main  memory  of  the  minicomputer  limited  the  total  number  of 
samples  taken  in  a  continuous  stream  to  18,432.  The  anemometer  signals  were  sampled 
at  3906  Hz  (6144  total  samples)  per  channel  or  2604  Hz  (9216  total  samples)  per 
channel  for  2  or  3  wire  applications,  respectively.  The  sampling  rates  resulted  in 
total  continuous  sample  periods  of  2.36  sec  for  both  2  and  3  wire  applications.  The 
digitized  voltage  samples  were  stored  on  magnetic  disks  using  a  DEC  Model  RX02  floppy 
disk  recorder. 

The  RX02  is  a  "double  density"  system  and  can  record  up  to  512  K  bytes  of  infor¬ 
mation  on  a  single  floppy  disk.  Reduction  of  the  voltage-time  records  to  either 
velocity-time  records  (isothermal  flow  -  2  wires)  or  velocity-temperature-time  records 
(flows  with  wall  heating  -  3  wires)  was  accomplished  off-line  using  an  LSI  11-03 
minicomputer.  The  reduced  temperature-velocity-time  or  velocity-time  results  were 
written  onto  double-density  magnetic  disks  and  copied  onto  magnetic  tape  for  purposes 
of  plotting  and  tabulation. 


4.  Data  Analysis  Techniques 

The  digitized  voltage  vs.  time  records  from  the  multi-wire  probes  were  reduced 
to  turbulence  quantities  using  an  LSI  11-03  minicomputer.  For  this  reduction  step 
the  digitized  data  were  read  into  the  computer  from  the  RX02  unit  while  the  following 
constants  for  each  sensor  were  input  through  a  terminal. 

R32  -  sensor  resistance  at  32°F 

Radiabatic  -  sensor  resistance  in  flow  with  no  overheat 
Rhot  ”  sensor  resistance  at  operating  temperature 

Rs  -  probe,  cable  and  anemometer  (40  n  for  TSI-1050  sets)  series  resistance 
a  -  temperature-resistance  coefficient 
A2»  B2  -  calibration  constants  from  Eq.  3 
k  -  calibration  constant  from  Eq.  4 

4.1  Isothermal  Flows  (2  Wire  Probes) 

Solution  for  the  velocity  components  (u  and  v  for  the  vertical  wire  arrays,  u 
and  w  for  the  horizontal  arrays)  for  each  time  step  proceeded  as  follows.  First, 
using  the  adiabatic  resistances  (no  sensor  overheat)  measured  for  the  sensors  in  the 
test  flow  and  Eq.  1,  the  temperature  (T)  of  the  flow  was  computed.  The  hot  sensor 
temperatures  (Tw)  were  then  computed  from  R^  and  Eq.  1.  For  each  time  step  the 
voltages  for  each  of  the  sensors  were  input  to  Eq.  3  to  determine  the  sensor  effec¬ 
tive  velocity  (Ue).  Next,  assuming  that  the  wires  of  the  x  array  were  perpendicular 
to  each  other  and  at  ±  45°  to  the  mainstream  flow  direction,  the  simultaneously 
measured  effective  sensor  velocities  were  combined  using  Eq.  4  to  find  either  u  and 
v  (vertical  array)  or  u  and  w  (horizontal  arrays).  As  a  check  on  the  accuracy  of 
the  assumption  that  the  wires  were  at  exactly  ±  45°  to  the  mainstream  one  of  the 
probes  was  also  calibrated  using  the  "wire  effective  angle"  method  of  Refs.  14  and  15. 


Voltage  vs.  time  records  were  reduced  to  fluctuating  velocity  components  using  these 
two  different  calibration-reduction  techniques  and  the  results  were  in  very  close 
agreement.  Once  the  velocity  component  vs.  time  record  was  generated  it  remained  a 
straightforward  matter  to  compute  any  desired  statistical  quantities  for  the  entire 
time  record.  The  following  turbulence  quantities  were  computed  for  the  u-v  (vertical 
array)  components.  Similar  quantities  with  the  transverse  velocity  component  (w) 
substituted  for  the  vertical  component  (v)  were  computed  for  the  horizontal  probe 
arrays. 


u,  u'2.  u7^,  u'4 

v,  v'2,  v'3,  v'4 
u'v' ,  u'2v' ,  u 1 v 1 2 
Su,  Sv,  Fu,  Fv 


-  first  through  the  fourth  moments  of  the  streamwise 
fluctuating  velocity 

-  first  through  the  fourth  moments  of  the  normal 
fluctuating  velocity 

-  double  and  triple  cross-products  (and  their  correlation 
coefficients 

-  skewness  and  flatness  of  both  velocity  components 


4.2  Flows  with  Mall  Heating  (3  Wire  Probes) 


For  the  3  wire  probes  employed  for  these  measurements  the  two  parallel  wires  of 
the  array  were  operated  at  different  overheats  (Rw-hot/Radiabatic  =  and  l-2)- 
The  data  reduction  technique  used  for  these  measurements  was  based  upon  the  assump¬ 
tion  that  during  any  time  step  the  effective  velocities  over  the  two  parallel  wires 
were  equal  (for  velocity  scales  equal  to  or  larger  than  the  separation  distance  be¬ 
tween  the  wires).  The  solution  technique  proceeded  by  first  finding  the  fluid  tem¬ 
perature  (T)  for  a  given  time  step.  Using  the  voltages  (Ew)  from  the  two  parallel 
sensors  and  assuming  that  Ue  was  equal  for  both  wires,  Eq.  3  was  iteratively  solved 
for  T.  Once  T  was  known  the  solution  for  the  velocity  components  (u  and  v)  for  each 
step  proceeded  as  in  4.1.  In  addition  to  computing  the  turbulence  quantities  listed 
in  section  4.1,  the  following  items  were  determined  for  the  cases  with  wall  heating. 

t,  t ' 2 ,  t'3#  t'4  -  first  through  the  fourth  moments  of  the  fluctuating 

_  _  temperature 

v't',  v't'2  -  velocity-temperature  cross-products  (and  their  corre¬ 

lation  coefficients) 

Sj,  Fy  -  skewness  and  flatness  of  temperature 


4.3  Reynolds  Stress  and  Turbulent  Heat  Flux  Corrections  for  Sensor  Separation 

The  accuracy  of  cross-products  of  correlated  turbulent  quantities  is  directly 
influenced  by  the  spacing  between  the  sensors  used  to  measure  these  quantities 
(Ref.  13).  The  contributions  of  the  smallest  scales  of  the  turbulence  (smaller  than 
the  transverse  sensor  spacing)  are  excluded  from  the  correlated  products.  As  examples 
of  the  impact  of  this  effect  Refs.  9,  10,  11,  12  and  14  all  present  Reynolds  stress 
measurements  (-u'v')  which  are  about  30  percent  lower  than  expected.  Unfortunately 
this  effect  cannot  be  eliminated  completely  because  a  minimum  limit  for  sensor  sepa¬ 
ration  is  reached  when  sensor  "cross-talk"  errors  become  significant.  As  discussed 
in  Section  2.1,  the  probes  used  for  the  present  program  were  specifically  designed 
to  minimize  these  effects. 
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Correction  factors  for  the  cross-product  terms  measured  in  this  program  were 
determined  by  the  following  technique.  Using  the  parallel  wires  of  the  3  wire  probes, 
the  transverse  spatial  correlation  coefficient  parallel)  was  determined  as  a  func¬ 
tion  of  position  (y/6)  in  the  test  boundary  layers.  Wire  separation  distances  (r) 
for  each  of  the  probes  were  accurately  measured  using  a  Nikon  Model  II  toolmaker's 
microscope.  By  assuming  that  the  correlation  coefficient  fell  with  the  square  of 
the  separation  distance  (r^)  (Ref.  16)  an  appropriate  spatial  correlation  coeffi¬ 
cient  could  then  be  calculated  for  any  x-type  probe/ boundary  layer  location  combina¬ 
tion.  Next,  assuming  that  the  contributions  to  the  cross-products  were  directly 
proportional  to  the  spatial  correlation  coefficient,  a  correction  factor  for  the 
probe/location  combination  was  determined: 


/J:  ('“^poroiieil 


r  (x  probe)]2 
r(poroltei)  I 


(5) 


As  an  example  the  correction  procedure  for  a  measured  Reynolds  stress  was  as  follows: 


i~U7T\  .  V  ^meosured 

'  u  v  'corrected  TTa 


(6) 


Typical  correction  factors  (b)  for  the  various  probes,  quantities  and  locations 
ranged  from  0.12  to  0.2.  A  journal  article  documenting  the  development  of  this  cor¬ 
rection  technique  is  currently  in  preparation. 


12 


EXPERIMENTAL  DATA 


1 .  Experimental  Test  Program 

Measurements  were  obtained  for  three  test  flow  conditions  of  incompressible, 

zero  pressure  gradient  flow  along  a  flat,  uniform  heat  flux,  test  wall.  The  three 

test  cases  of  this  program  reproduced  conditions  employed  for  an  earlier  AFOSR  Con¬ 
tract  (Ref.  2).  For  all  test  cases  the  free-stream  velocity  was  100  ft/sec  and  the 
test  surface  boundary  layer  passed  through  natural  transition,  i.e.,  no  artificial 
trips  were  employed  to  promote  boundary  layer  transition.  Data  were  obtained  for 
three  levels  of  free-stream  turbulence,  (1)  at  the  tunnel  minimum  turbulence  level 
and  (2)  and  (3)  at  higher  levels  of  free-stream  turbulence  generated  with  bi-plane 
grids.  Using  the  nomenclature  of  Refs.  1  and  2,  the  free-stream  turbulence  levels 
of  this  program  are  designated  as  follows 

(1)  No  grid-low  free-stream  turbulence  (Grid  0)  -  Te  s  \  percent 

(2)  Grid  number  2  (%  in.  bars)  -  Tenom  =  2  percent 

(3)  Grid  number  4  (2  in.  bars)  -  Tenom  =  6  percent 

A  complete  documentation  of  the  multi-component  turbulence  decay,  integral 
length  scale  growth  and  spectral  distributions  generated  by  these  particular  test 
grids  is  available  in  Refs.  1  and  2. 

For  each  of  these  flow  conditions  experimental  boundary  layer  profile  data  were 
obtained  at  three  streamwise  locations  (x  s  52,  68  and  84  inches)  for  both  an  adia¬ 
batic  test  surface  (no  wall  heating)  and  with  a  uniform  surface  heat  flux  condition. 
With  no  wall  heating  the  following  data  were  measured: 


Type  of  Data 

Instrumentation 

Measurement  Stations 
Per  Profile 

Profiles  of  streamwise 
velocity 

(mean  and  fluctuating) 

Single,  horizontal,  linearized 
hot  wire 

30 

Profiles  of  streamwise  and 
normal  velocities 

Vertical  x-type  wires  with 
analog-digital  data  system 

17 

Profiles  of  streamwise  and 
transverse  velocities 
(mean  and  fluctuating) 

Horizontal  x-type  wires  with 
analog-digital  data  system 

17 

With  the  uniform  wall  heat  flux  conditions  the  following  data  were  measured: 
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Type  of  Data 

Instrumentation 

Measurement  Stations 
Per  Profile 

Surface  Stanton  Number 
distribution 

Thermocouple  instrumentation 
incorporated  into  uniform  heat 
flux  test  surface 

210 

surface  locations 

Profiles  of  streamwise 
mean  velocity 

Miniature  boundary  layer  pitot 
probes 

90 

Profiles  of  mean 
temperature 

Miniature  boundary  layer  ther¬ 
mocouple  probes 

90 

Profiles  of  temperatures 
and  streamwise  and  normal 
velocities 

(mean  and  fluctuating) 

3-Wire  probes  with  analog- 
digital  data  system 

— 

17 

In  summary,  for  each  of  che  three  flow  conditions  surface  heat  transfer  distri¬ 
butions  and  three  stations  of  profile  data  were  measured.  In  total  (3  conditions)  x 
(3  profiles)  x  (6  types  of  profile  data)  =  54  profile  surveys  were  documented. 


2.  Boundary  Layer  Profile  Data  Format 

The  mean  and  fluctuating  quantities  measured  for  the  various  flow  conditions  and 
profile  locations  have  been  assembled  in  both  graphical  and  tabular  form.  Comparisons 
of  tnese  results  for  the  various  flow  conditions  and  with  similar  results  from  other 
experiments  will  be  presented  in  the  Analysis  of  Results  section  below.  The  compiled 
data  for  all  the  measurements  stations  are  given  in  Appendix  B  -  Experimental  Data. 

As  a  guide  to  the  format  of  the  presentation  and  results,  the  data  for  a  single  sample 
profile  are  given  in  Figs.  4  through  5E  and  in  Tables  1  through  3B.  These  particular 
sample  profile  data  were  obtained  at  the  "middle"  free-stream  turbulence  level  (Grid 
2,  Te  =  1.6%)  at  X  =  68  in.  The  mean  profile  (total  pressure  and  thermocouple  probe) 
data  for  the  sample  set  are  presented  in  graphical  form  in  Fig.  4  and  in  tabular  form 
in  Table  1.  Table  2  presents  a  compilation  of  the  test  flow  conditions  and  values 
computed  from  the  mean  profile  data.  These  mean  profile  data  are  presented  both  in 
the  form  of  velocity  and  temperature  ratios  versus  y/6  and  in  the  coordinates  of  the 
universal  velocity  and  temperature  "laws  of  the  wall".  As  discussed  in  the  Data 
Analysis  Techniques  section,  the  digitized  data  reduction  system  made  possible  the 
computation  of  any  desired  moments  and  cross-products  of  the  various  measured  fluc¬ 
tuating  quantities.  A  total  of  38  quantities  were  selected  for  presentation  on  the 
grounds  that  they  met  either  or  both  of  the  following  criterion:  (1)  the  quantity 
is  employed  in  some  existing  boundary  layer  turbulence  modeling  method  or  (2)  the 
quantity  serves  as  a  diagnostic  of  the  characteristics  of  the  turbulence,  e.g., 
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intermittency.  These  various  quantities  were  both  plotted  and  tabulated  for  each 
profile  station  (see  Appendix  8).  The  results  for  the  sample  profile  are  presented 
in  Figs.  5A  through  E  and  in  Table  3A  and  B.  The  distribution  of  a  series  of  turbu¬ 
lence  quantities  computed  from  the  fluctuating  velocity  data  are  presented  in  Fig. 
5A.  Starting  in  the  upper  left-hand  corner  of  the  figure,  the  distributions  of  the 
individual  components  of  the  turbulence  are  compared  with  the  results  of  Klebanoff 
(obtained  for  near-zero  free-stream  turbulence).  Moving  clockwise,  the  next  figure 
presents  the  measured  distributions  of  the  Reynolds  stress  nondimensional ized  by  the 
friction  velocity  (UT  was  determined  independently  from  the  mean  profile  data)  and 
the  distribution  of  the  shear  stress  correlation  coefficient.  Also  included  in  this 
plot  is  the  distribution  of  shear  stress  computed  from  the  mean  profile  data  using 
the  technique  of  Ref.  17.  The  lower  right-hand  corner  plot  of  Fig.  5A  presents  the 
transport  velocities  of  turbulent  shear  stress  and  kinetic  energy  defined  as  follows 


uV‘ 

uV 
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See  Ref.  14,  pp.  220-239,  and  Ref.  18  for  the  development  of  these  terms.  (Due  to  a 
software  error  VT  was  not  computed  for  the  data  of  Grid  2  and  does  not  appear  for 
this  sample  plot.)  The  remaining  plot  of  Fig.  5A  presents  the  structural  coefficient 
(as  defined  in  Ref.  19)  distributions  for  this  case. 

Turbulence  quantities  computed  from  the  fluctuating  velocities  and  temperatures 
are  presented  in  Fig.  5B.  Distributions  of  the  turbulent  heat  flux  and  its  correla¬ 
tion  coefficient  are  given  for  the  plot  in  the  upper  left-hand  corner.  The  turbulent 
heat  flux  distribution  is  shown  nondimensional ized  by  the  independently  measured  wall 
heat  flux.  Also  shown  in  this  plot  is  the  distribution  of  heat  flux  through  the  boun¬ 
dary  layer  as  computed  from  the  mean  velocity  and  temperature  profile  data  (see  Ref. 
17).  The  upper  right-hand  corner  plot  of  Fig.  5B  presents  the  fluctuating  tempera¬ 
ture  distribuitons  in  two  forms:  (1)  nondimensional ized  by  the  friction  temperature 
and  (2)  nondimensional ized  by  the  temperature  difference  between  the  wall  and  free- 
stream.  The  lower  two  plots  of  Fig.  5B  give  the  distributions  of  the  turbulent 
Prandtl  number,  Prt,  and  two  structural  coefficients,  a-|^e  and  G-j  e  (see  Ref.  19). 
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Triple  product  distributions  of  the  fluctuating  velocity  components  are  given 
in  Fig.  5C.  Streamwise-transverse  (u'w*)  products  are  grouped  on  the  left-hand  plot 
while  streamwise-normal  (u'v')  products  appear  in  the  right. 


For  Fig.  5D  the  left-hand  figure  presents  the  skewness  factor  distributions  for 
the  fluctuating  temperatures  and  velocities.  Note  that  skewness  factor  distributions 
of  the  streamwise  component  (u‘)  were  determined  both  from  the  data  from  the  vertical 
x  probes  (Suv)  and  the  horizontal  x  probes  (Su^).  The  correlation  coefficients  for 
the  triple  products  of  Fig.  5C  are  given  in  the  right-hand  plot  of  Fig.  50. 


Flatness  factor  distributions  for  the  fluctuating  temperatures  and  velocities 
are  given  in  Fig.  5E.  As  with  the  skewness  factors  of  Fig.  5D  flatness  factor  dis¬ 
tributions  of  the  streamwise  component  are  given  for  both  the  data  for  the  vertical 
x  (FUv)  and  horizontal  x  (F^)  probes.  To  avoid  crowding  on  the  figure  the  flatness 
of  the  temperature  fluctuations  was  plotted  after  dividing  by  2. 


Tabulated  values  of  these  fluctuating  quantities  are  given  in  Table  3A  and  3B. 
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ANALYSIS  AND  DISCUSSION  OF  EXPERIMENTAL  RESULTS 


The  present  experimental  test  program  was  designed  to  examine,  in  detail,  the 
effects  of  free-stream  turbulence  on  the  heat  transfer  through  turbulent  boundary 
layers.  The  test  conditions  for  the  present  program  were  intended  to  reproduce  cases 
for  which  other  experimental  data  had  been  obtained  previously  under  an  earlier 
AFOSR  contract  (Ref.  2). 


1.  Comparisons  of  Present  Results  with  Results  of  Ref.  2 

A  number  of  the  measurements  reported  in  Ref.  2  were  repeated  during  the  present 
program  providing  a  measure  of  consistency  for  the  two  sets  of  data.  Note  that  these 
various  measurements  were  obtained  in  the  same  test  facility  and  on  the  same  test 
surface  but  about  three  years  apart. 

1.1  Free-Stream  Turbulence  Data 


Measurements  of  the  components  of  the  free-stream  turbulence  (outside  the  boun¬ 
dary  layer)  were  generally  in  excellent  agreement  (absolute  levels  of  individual 
components  agreed  within  0.3  percent)  with  the  measurements  of  Ref.  2.  The  excep¬ 
tions  to  this  rule  were  the  transverse  (w‘)  fluctuation  levels  measured  for  Grid  4. 
These  data  were  consistently  higher  (relatively  15%  higher)  than  those  measured  for 
the  same  flow  condition  of  Ref.  2.  This  discrepancy  will  be  discussed  in  more  detail 
in  a  following  section. 

1.2  Heat  Transfer  Distributions 

For  all  three  test  conditions  the  agreement  between  the  Stanton  numbers  measured 
for  the  present  program  and  for  Ref.  2  was  excellent  (±1%).  The  heat  transfer  dis¬ 
tributions  for  the  no  grid  (Grid  0),  Grid  2  and  Grid  4  cases  were  virtually  identical 
with  those  presented  in  Ref.  2  in  Figs.  41,  49  and  57,  respectively  and  are  not  re¬ 
peated  here. 

1.3  Boundary  Layer  Transition  Location 

For  the  Grid  2  and  4  test  cases  the  agreement  between  the  present  transition 
location  data  and  the  similar  data  of  Ref.  2  was  within  ±  3  percent.  For  the  no-grid 
case,  however,  the  transition  Reynolds  number  increased  from  Rex  =  1.2  x  106  to  Rex  = 
1.35  x  106.  This  change  in  observed  transition  location  was  related  to  the  three 
dimensional  character  of  the  transition  process  for  the  no-grid  case.  As  discussed 
in  Refs.  1  and  2,  test  section  corner  flows  contaminate  the  flat  test  wall  laminar 
boundary  layer  for  the  low  free-stream  turbulence  case  and  produce  premature  transi¬ 
tion  along  the  tunnel  centerline.  This  sidewall  contamination  was  not  important 
for  the  higher  levels  of  free-stream  turbulence  because  two-dimensional  natural 
transition  resulted  well  upstream  of  these  effects.  For  the  present  no-grid  test 
conditions  the  leading  edge  scoop  adjustment  was  improved  over  the  setting  of  the 


tests  of  Ref.  2  producing  reduced  secondary  corner  flows  and  an  increased  transition 
Reynolds  number.  The  transition  Reynolds  number  (Ree)  for  the  present  tests  was  in 
excellent  agreement  with  classic  two-dimensional  transition  vs.  turbulence  correla¬ 
tions.  The  turbulent  boundary  layers  which  developed  downstream  of  transition  for 
both  the  present  test  and  the  test  of  Ref.  2  both  exhibited  the  classic  characteris¬ 
tics  for  zero-pressure  gradient,  low  free-stream  turbulence,  two-dimensional  flow. 

1 .4  Mean  Profile  Data 

Integral  thicknesses  (<5*  and  e)  computed  from  the  prifle  data  for  the  Grid  2  and 
4  cases  agreed  within  2  percent  with  the  results  from  Ref.  2  at  the  respective 
locations.  For  the  no-grid  case  the  integral  thicknesses  were  about  12  percent  re¬ 
duced  from  those  computed  for  Ref.  2  at  similar  stations.  This  reduction  in  boundary 
layer  thickness  resulted  from  the  increased  length  of  laminar  flow  upstream  of  tran¬ 
sition  for  the  present  no-grid  data. 

If  comparisons  are  made  only  for  profiles  with  equal  Ree  the  results  for  all  the 
grids  (0,  2  and  4)  are  practically  identical  to  the  respective  cases  in  Ref.  2.  When 
plotted  in  U+  or  T+  vs.  Y+  coordinates  the  mean  velocity  and  temperature  profiles  ex¬ 
hibit  significant  regions  (30  <  Y+  <  300)  of  logarithmic  behavior.  Both  the  velocity 
and  temperature  wakes  showed  significant  reduction  with  increasing  free-stream  turbu¬ 
lence  as  did  the  similar  data  of  Ref.  2.  Skin  friction  coefficients  computed  from 
fits  of  the  mean  velocity  data  to  the  "law-of-the-wall "  were  in  excellent  agreement 
with  the  results  of  Ref.  2. 

1.5  Comparison  with  Earlier  Data  -  Conclusion 

The  conclusion  reached  from  the  free-stream  turbulence,  heat  transfer  and  mean 
profile  data  was  that  the  test  conditions  of  Ref.  2  were  closely  duplicated  for  the 
present  series  of  measurements.  In  effect  these  present  measurements  can  be  con¬ 
sidered  as  an  additional  set  of  data  for  the  same  test  conditions  as  Ref.  2. 


2.  Profile  Data  with  Low  Free-Stream  Turbulence 

A  number  of  comparisons  have  been  made  between  the  data  obtained  for  the  present 
j  no-grid  (low  free-stream  turbulence)  profiles  and  measurements  from  other  experiments. 

I  These  comparisons  are  intended  to  provide  a  measure  of  the  accuracy  and  consistency 

of  the  present  boundary  layer  turbulence  data. 

*  The  distributions  of  the  u'  and  v1  components  of  the  turbulence  profiles  were 

i  in  very  good  agreement  with  the  classic  results  of  Klebanoff  (Ref.  20),  see  for 

example  Appendix  B— Fig.  B-4A.  The  transverse  component  (w’)  measurements,  however, 
were  typically  about  15  percent  reduced  from  Klebanoff' s  results  with  w1  only 
slightly  greater  than  v'.  These  present  w'  distributions  are  thought  to  be  accurate 
as  they  are  in  close  agreement  (as  were  the  u'  and  v'  distributions)  with  the  recent 
results  of  Ref.  14.  The  turbulent  shear  stress  distributions  measured  for  these 
low  free-stream  turbulence  cases  were  in  excellent  agreement  with  the  shear  stress 


distributions  computed  from  the  mean  profiles.  Fig.  B-4A.  The  accuracy  of  these 
particular  measurements  is  also  supported  by  the  fact  that  for  all  cases  the  u'v' 
correlation  coefficient  was  near  the  widely  accepted  value  of  0.44  (Ref.  21)  across 
most  of  the  boundary  layer.  The  measured  distributions  of  the  "structural "  coeffi¬ 
cients  were  in  good  agreement  with  the  widely  accepted  constants,  ai  =  u'v'/q2  = 

0.15,  a2  =  ITW-  =  0.5,  a3  =  =  0.2. 

A  number  of  the  turbulence  quantities  computed  from  the  present  data  were  also 
determined  for  a  similar  zero-pressure  gradient,  low  free-stream  turbulence  flow  in 
the  work  of  Ref.  14.  The  present  distributions  of  transport  velocity  of  turbulent 
shear  stress  (VT),  transport  velocity  of  turbulent  kinetic  energy  (Vg)  (Fig.  B-4A), 
u'v'  and  u'w'  triple  products  (Fig.  B-4C),  skewness  factors  (Fig.  B-4D)  and  flatness 
factors  (Fig.  B-4E)  were  all  in  good  agreement  with  the  respective  data  of  Ref.  14. 

The  turbulent  heat  flux  distribution  measurements  were  in  reasonably  close 
agreement  with  the  distributions  inferred  from  the  mean  profiles  (Fig.  B-4B).  The 
accuracy  of  these  mean  profile  distributions  is  unclear  because  of  extreme  sensi¬ 
tivity  to  uncertainties  in  the  mean  temperature  profiles.  The  fluctuating  tempera¬ 
ture  distributions  agreed  very  closely  with  the  distributions  measured  in  Ref.  21 
(Fig.  B-4B).  The  values  of  the  thermal  coefficient  a]t  were  about  30  percent  greater 
than  those  determined  in  Ref.  21  with  the  cause  of  the  difference  uncertain.  The 
authors  were  unaware  of  any  other  measurements  of  G-j0  to  which  the  present  data 
could  be  compared.  Finally  the  turbulent  Prandtl  number  distributions  measured  for 
the  low  free-stream  turbulence  cases  were  in  excellent  agreement  with  the  proposed 
distribution  of  Rotta  (Ref.  22). 

There  were,  then,  a  large  number  of  experimentally  determined  turbulence  quanti¬ 
ties  in  the  present  program  which  agreed  very  closely  with  the  results  of  other 
studies.  The  conclusions  reached  from  this  result  are  that  one  can  have  a  high  level 
of  confidence  in  the  present  data  acquisition  and  reduction  technique  and  that  the 
turbulence  quantities  reported  here  can  be  expected  to  be  both  accurate  and  self- 
consistent. 


3.  Effects  of  High  Free-Stream  Turbulence 
on  the  Fluctuating  Velocities 

The  impact  of  increased  free-stream  turbulence  on  the  boundary  layer  turbulence 
kinetic  energy  distribution  is  shown  in  Fig.  6.  Presented  in  this  figure  are  ex¬ 
perimental  data  from  the  present  program  obtained  at  stations  with  nearly  equal 
momentum  thickness  Reynolds  numbers  (Re0  *  5500  ±  100)  for  Grids  0,  2  and  4.  Also 
given  in  Fig.  6  are  turbulence  kinetic  energy  distributions  measured  for  similar 
zero  pressure  gradient,  low  speed  flows  by  other  investigators  (Refs.  12  and  21). 
Integral  thicknesses  were  not  computed  for  these  other  data  but  based  upon  the 
stated  values  of  6  it  is  estimated  that  for  Ref.  12  Re©  was  about  3500  while  for 
Ref.  21  Re0  ~  700.  Agreement  between  the  present  low  free-stream  turbulence  profile 
case  (Grid  0)  and  the  similar  data  of  Ref.  21  is  reasonably  good  except  very  close 
to  the  wall.  As  discussed  earlier  this  near-wall  discrepancy  resulted  from  the 
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relatively  higher  values  of  w'  determined  in  the  study  of  Ref.  21.  For  the  profiles 
with  higher  free-stream  turbulence  there  were  no  cases  where  the  present  flow  condi¬ 
tions  (Re^  and  Te)  were  identical  with  those  of  Ref.  12.  Despite  this  a  comparison 
of  these  results  shows  reasonably  good  agreement  for  both  the  trends  and  abs^1  +:e 
magnitudes  of  the  kinetic  energy  distributions.  Both  the  present  data  and  the  inde¬ 
pendent  results  from  Refs.  12  and  21  indicate  a  progressive  increase  in  boundary 
layer  turbulence  kinetic  energy  with  increasing  free-stream  turbulence.  Increased 
levels  of  turbulence  kinetic  energy  were  measured  across  the  entire  thickness  of  the 
boundary  layer. 

The  effects  of  the  free-stream  turbulence  level  on  the  indivdiual  components  of 
the  boundary  layer  turbulence  for  these  same  three  profiles  are  shown  in  Fig.  7.  The 
streamwise  (u1)  component  followed  the  same  trends  as  the  turbulence  kinetic  energy, 
increasing  with  free-stream  turbulence  level  over  the  entire  thickness  of  the  boun¬ 
dary  layer.  The  normal  component  (v‘),  however,  was  damped  by  the  presence  of  the 
solid  wall  and  showed  virtually  no  change  over  the  lower  half  of  the  boundary  layer. 
See  Ref.  23  for  an  in-depth  study  of  the  interaction  of  solid  surfaces  with  turbulent 
fluctuations.  The  distribution  of  w1  for  Te  =  4.2  percent  showed  a  large  increase 
over  the  distributions  for  the  lower  turbulence  levels.  Outside  the  boundary  layer 
w1  was  also  measured  to  be  about  20  percent  higher  than  u'  or  v'.  On  the  grounds 
that  earlier  independent  measurements  of  the  free-stream  turbulence  for  this  grid 
showed  the  turbulence  to  be  isotropic  at  this  station  (Ref.  1)  and  the  unreasonably 
large  "jump"  in  the  w'  distribution  across  the  entire  boundary  layer  it  has  been 
concluded  that  these  w’  measurements  are  in  error.  It  is  thought  that  there  was  an 
error  in  the  calibration  for  the  horizontal  x  wire  probe  used  for  the  Grid  4  test 
cases.  The  Grid  4  w'  data  are  reported  here  as  measured,  that  is  uncorrected  for 
this  probable  error.  It  is  estimated  that  by  reducing  the  measured  w'  data  by  20 
percent  a  reasonably  accurate  set  of  distributions  of  the  transverse  component  for 
this  Grid  4  case  would  result. 

With  the  Grid  4  w1  data  reduced  by  20  percent  the  conclusion  that  can  be  reached 
from  Fig.  7  is  that  the  u‘  and  w‘  component  of  turbulence  increased  progressively  with 
increasing  free-stream  turbulence  level.  Both  components  increased  at  all  locations 
in  the  boundary  layer.  The  vertical  component  v',  however,  was  essentially  constant 
and  independent  of  free-stream  turbulence  level  for  the  lower  half  of  the  boundary 
layer. 

Distributions  of  the  boundary  layer  turbulence  structural  coefficients  (Bradshaw, 
et  al..  Ref.  19)  are  given  in  Figs.  8A  and  8B.  As  can  be  seen  from  an  inspection  of 
Fig.  8A  the  ratio  of  shear  stress  to  turbulence  kinetic  energy  (a-j)  decreased  across 
the  entire  boundary  layer  with  increasing  free-stream  turbulence.  The  observed 
decrease  was  most  extreme  over  the  outer  60  percent  of  the  boundary  layer.  Also 
shown  with  the  present  a-|  distribtuion  data  are  similar  results  from  Refs.  14  and  21. 
Agreement  between  these  similar  (not  identical,  Ree  and  Te  were  slightly  different) 
sets  of  data  was  very  good. 

The  ratios  of  the  direct  stress  components  to  the  turbulent  kinetic  energy 
(a2»  a3  and  84)  are  given  in  the  remaining  plots  of  Figs.  8A  and  8B.  Employing  the 
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previously  described  20  percent  reduction  to  the  w‘  component  for  Te  =  4.2  percent 
(this  also  reduces  q?  for  Te  =  4.2  percent),  fairings  of  the  corrected  structural 
coefficient  distributions  for  the  highest  turbulence  level  are  given  in  the  figures. 
Using  the  measured  results  for  Te  *  0.2  percent  and  1.5  percent  and  the  corrected 
fairings  for  Te  =  4.2  percent  the  following  conclusions  were  reached.  As  the  free- 
stream  turbulence  stress  level  was  increased  a2  (u'Z/q?)  increased  slightly  above 
0.5,  the  value  widely  used  for  low  free-stream  turbulence  boundary  layers.  The 
greatest  percentage  change  was  observed  for  the  lower  half  of  the  boundary  layer  for  a3 
(a3  =  v'2/q2).  This  ratio  decreased  progressively  with  increasing  Te  dropping  to 
about  0.12  (40%  reduction  from  classic  value  of  0.2)  for  Te  =  4.2  percent.  Only 
very  small  changes  were  observed  for  54  (a4  =  w'2/qZ)  with  the  measured  values 
grouping  around  0.3  for  the  lower  half  of  the  boundary  layer. 

A  number  of  previous  studies  of  free-stream  turbulence  effects  on  turbulent 
boundary  layers  (Refs.  8,  10,  11,  14  and  234)  have  reported  finite  turbulent  shear 
stress  levels  beyond  the  edge  of  the  velocity  boundary  layer.  This  effect  was  also 
observed  for  the  present  program.  A  comparison  between  the  present  results  and 
those  of  the  previous  investigations  is  given  in  Fig.  9  where  the  turbulent  shear 
level  at  the  edge  of  the  boundary  layer  (6  0.995)  is  given  as  a  function  of  Te.  The 
data  from  the  present  study  and  the  results  reported  for  most  of  the  other  experiments 
are  tightly  grouped.  Taken  together,  these  data  indicate  an  increase  of  turbulent 
shear  at  the  boundary  edge  directly  proportional  to  the  free-stream  turbulence  level. 
Huffman's  results,  which  are  believed  (Ref.  2)  to  contain  significant  errors  due  to 
anistropy,  show  much  larger  levels  of  turbulent  shear  than  the  other  studies. 

Additional  evidence  of  the  impact  of  free-stream  turbulence  on  the  characteris¬ 
tics  of  the  turbulence  near  the  boundary  layer  edge  is  provided  by  the  measurements 
of  flatness  factor.  The  flatness  factor  (u^/(u'2)2  is  an  indication  of  the  dis¬ 
tribution  of  velocity  fluctuations  in  a  set  of  samples.  For  a  normal  Gaussian  dis¬ 
tribution  the  flatness  factor  is  equal  to  3  with  larger  values  indicating  contri¬ 
butions  from  intermittent  turbulent  fluctuations.  Flatness  factor  distributions  of 
the  streamwise  fluctuating  velocity  component  are  given  in  Fig.  10  as  a  function  of 
position  in  the  boundary  layer.  An  examination  of  Fig.  10  shows  that  the  inter¬ 
mittent  character  of  the  turbulence  near  the  edge  of  the  boundary  layer  was  greatly 
reduced  by  increased  free-stream  turbulence.  For  a  boundary  layer  beneath  a  low 
turbulence  mainstream  a  relatively  sharp  irregular  "edge"  of  turbulent  boundary 
layer  flow  results  adjacent  to  the  non-turbulent  freestream.  With  higher  levels  of 
free-stream  turbulence  this  distinct  border  appears  to  have  disappeared. 


4.  Effects  of  High  Free-Stream  Turbulence 
on  the  Turbulent  Prandtl  Number 

The  measured  distributions  of  turbulent  shear  stress  (u'v'),  the  turbulent  heat 
flux  (v't1)  and  the  normal  derivatives  of  the  mean  velocity  and  temperature  were 
combined  to  form  local  turbulent  Prandtl  numbers. 


Pr 


♦ 


«h 


21 


v  — 


02) 


where  em  =  eddy  diffusivity  of  mass 
eh  =  eddy  diffusivity  of  heat 


_  dT 

-uV  dy 

v7?” 

3y 


(9) 


The  distributions  of  turbulent  Prandtl  number  measured  for  the  various  test 
cases  are  presented  in  Fig.  11.  The  results  from  all  three  profile  locations  (x  = 
52,  68  and  84  inches)  for  all  three  free-stream  turbulence  levels  are  included  in 
Fig.  11  with  an  average  free-stream  turbulence  level  assigned  to  each  set.  For  all 
points  above  the  wall  the  turbulent  Prandtl  numbers  were  determined  from  the  turbu¬ 
lent  heat  flux  and  shear  stresses  measured  with  the  hot  wire  probes  and  from  the 
derivatives  of  the  mean  profiles  measured  with  the  total  pressure  and  thermocouple 
probes.  At  the  wall  the  turbulent  Prandtl  numbers  were  determined  from  the  mean 
temperature  and  velocity  profile  data  by  assuming  that  for  at  least  some  small  dis¬ 
tance  the  ratio  of  shear  stress  to  heat  flux  remains  at  the  wall  value. 
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Near-wall  values  of  the  turbulent  Prandtl  number  were  evaluated  from  Eq.  15 
using  friction  velocities  (Ux)  determined  from  the  mean  velocity  profile  fits  to  the 
"law-of-the-wall".  Values  of  3T/3U  were  determined  graphically  from  the  near-wall 
velocity  and  temperature  profile  data. 

Errors  in  the  four  measured  terms  of  Eq.  9  combined  to  produce  considerable 
scatter  in  the  data  of  Fig.  11.  This  scatter,  however,  is  much  less  than  reported 
for  the  similar  measurements  of  Refs.  21  and  25.  It  is  expected  that  the  consistency 
and  absolute  accuracy  of  such  local  turbulent  Prandtl  number  measurements  could  be 
further  improved  by  employing  larger  samples  of  the  turbulent  data. 

The  turbulent  Prandtl  number  distributions  measured  for  the  low  free-stream 
turbulence  profiles  were  in  good  agreement  with  the  similar  data  of  Ref.  17.  (Ref. 

17  employed  mean  profile  data  only.)  In  addition,  Rotta's  (Ref.  22)  suggested  Prt 
distribution  for  low  free-stream  turbulence  boundary  layers  appears  to  represent 
the  present  low  turbulence  data  well. 
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Pr,  =  0.95-0.45  (y/8)2 

An  examination  of  Fig.  11  indicates  that  as  the  free-stream  turbulence  level  was 
raised  the  turbulent  Prandtl  number  increased  over  nearly  the  entire  boundary  layer. 
Values  of  Pr^  over  unity  were  recorded  at  y/6  -  0.3  for  the  highest  free-stream  tur¬ 
bulence  level.  These  increased  outer  region  turbulent  Prandtl  numbers  for  high 
free-stream  turbulence  levels  had  not  been  expected.  At  the  outset  of  the  program 
the  turbulence  characteristics  of  the  outer  portion  of  the  boundary  layer  were  known 
to  be  altered  considerably  by  increased  levels  of  free-stream  turbulence.  It  was 
also  known  from  the  measurements  of  Ref.  1  that  the  free-stream  turbulence  had  a 
large  impact  on  the  turbulent  heat  transfer  with  the  Reynolds  analogy  factor  increas¬ 
ing  with  increasing  turbulence  level.  It  was  speculated  that  the  increased  Reynolds 
analogy  factor  might  result  from  lowered  Prt  levels  (relatively  greater  increase  in 
v't*  as  comapred  to  u'v')  in  the  outer  portions  of  the  boundary  layer.  The  experi¬ 
mental  results  of  Fig.  11  indicate  just  the  opposite  effect.  As  the  free-stream  tur¬ 
bulence  level  was  increased  the  outer  region  Prt  levels  increased  while  the  near- 
wall  values  (determined  from  the  mean  profile  data)  indicate  a  small  but  progressive 
decrease.  The  following  expression,  a  modification  of  Rotta's  (Eq.  16)  low  free-stream 
turbulence  equation,  represents  the  measured  results  reasonably  well. 


Pr,  =  (jo.95  -  0.45  (y/S)2]  (l  +  2T)2]- 
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Equation  17  is  shown  in  Fig.  11  for  the  three  turbulence  levels  for  which  the 
experimental  data  were  obtained.  At  Te  =  0.2  percent,  Eq.  17  is  practically 
identical  to  Rotta's  (Eq.  16)  expression. 

The  mean  velocity  and  temperature  profile  data  from  both  the  present  program  and 
from  Ref.  1  provide  additional  evidence  that  the  near-wall  turbulent  Prandtl  number 
decreased  with  increased  free-stream  turbulence  level.  (The  arguments  for  this  con¬ 
clusion  will  be  presented  here  in  a  highly  abbreviated  form.  A  more  in-depth  exami¬ 
nation  of  these  effects  will  be  conducted  during  the  preparation  of  a  technical 
journal  article  on  this  contract  work.)  The  effects  of  the  free-stream  turbulence 
on  the  similarity  between  the  mean  velocity  and  temperature  profiles  was  examined  by 
plotting  the  velocity  ratio  (U/Ue)  versus  the  temperature  ratio  (Tw-T/Tw-Te)  across 
the  boundary  layers.  Plots  of  these  ratios  for  all  profile  locations  and  free-stream 
turbulence  levels  are  given  in  Fig.  12.  Also  given  in  Fig.  12  are  the  similar  data 
from  the  same  stations  and  turbulence  levels  obtained  in  Ref.  1.  An  examination 
of  Fig.  12  indicates  that  for  all  cases,  independent  of  the  free-stream  turbulence 
level,  the  mean  velocity  and  temperature  profiles  remained  highly  similar.  This 
similarity  between  the  velocity  and  temperature  profiles  extends  across  at  least  the 
outer  90  percent  of  the  boundary  layer  thickness  including  all  the  wake  and  at  least 
some  of  the  logarithmic  zone.  It  follows  that  the  shapes  of  the  velocity  and  tem¬ 
perature  profiles  should  also  be  similar  when  plotted  in  universal  (U+  or  T+  vs.  Y+) 
coordinates.  It  has  been  observed  in  virtually  every  study  of  free-stream  turbulence 
effects  on  turbulent  boundary  layers  that  the  wake  strength  of  the  velocity  boundary 
layer  was  progressively  reduced  with  increasing  free-stream  turbulence.  It  was  also 
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observed  in  Refs.  2  and  26  that  for  a  given  turbulence  level  the  temperature  wake 
was  reduced  by  a  larger  amount  than  was  the  velocity  wake.  Implicit  in  the  formula¬ 
tion  of  the  temperature  law-of-the-wall  is  the  assumption  that  the  turbulent  Prandtl 
number  is  constant  across  the  entire  boundary  layer.  The  large  thermal  wake  depres¬ 
sions  reported  in  Refs.  2  and  26  followed  from  the  use  of  an  average  boundary  layer 
turbulent  Prandtl  number  for  all  the  profiles.  The  following  interpretation,  how¬ 
ever,  is  more  consistent  with  the  conclusion  from  Fig.  12,  that  the  shapes  of  the 
outer  region  velocity  and  temperature  profiles  remained  similar  for  all  turbulence 
levels  and  streamwise  locations.  If  the  near-wall  turbulent  Prandtl  numbers  were 
assumed  to  be  reduced  with  increased  Te  (as  Fig.  11  indicates)  the  slope  of  the  tem¬ 
perature  law-of-the-wall  ( 1 /k0  =  Prt/K)  would  be  reduced.  With  a  reduced  logarithmic 
region  slope  the  apparent  temperature  wake  strength,  which  is  the  maximum  deviation 
from  the  log-law,  would  increase.  An  examination  of  the  temperature  profiles  of 
Ref.  2  (in  T+  vs.  Y+  coordinates)  indicated  that  good  fits  to  the  temperature  law- 
of-the-wall  could  be  achieved  from  Y+  =  30  to  Y/6  =  0.1  if  nr,  was  set  equal  to  n* 

For  the  present  data,  very  good  agreement  between  the  thermal  and  velocity  wake 
strengths  resulted  from  the  use  of  the  near-wall  turbulent  Prandtl  numbers  of  Fig. 

11  for  the  respective  profiles. 

Finally,  with  regard  to  a  potential  physical  mechanism  producing  the  reduced 
near-wall  Pr^,  the  diffusion  terms  of  the  turbulence  kinetic  energy  transport  equa¬ 
tion  for  velocity  and  temperature  (Ref.  19)  differ  by  the  contribution  of  the  pressure- 
velocity  fluctuation  product.  Blom  (Ref.  27)  has  pointed  out  that  the  pressure  fluc¬ 
tuations  serve  to  transfer  energy  from  the  relatively  higher  u'  component  of  turbu¬ 
lence  to  the  relatively  smaller  v1  and  w'  components.  Blom  also  argued  that  the 
absence  of  the  p'v'  term  in  the  temperature  diffusion  term  could  explain  the  reduc¬ 
tion  of  Prj-  below  unity.  Since  the  effect  of  increased  Te  was  to  increase  the  diffe¬ 
rence  between  the  u'  and  v'  components  of  turbulence  near  the  wall  the  importance 
of  the  p'v'  term  may  grow  with  TP.  In  other  words,  with  increasing  difference  between 
u'  and  v'  the  effect  of  the  p'v'  term  may  be  to  progressively  decrease  Prj-. 

The  overall  impact  of  free-stream  turbulence  on  boundary  layer  heat  transfer 
rates  is,  then,  to  depress  the  near-wall  turbulent  Prandtl  number.  Since  the  heat 
and  momentum  transport  in  turbulent  boundary  layers  are  dominated  by  the  turbulent 
eddy  contributions  the  result  is  that  the  Reynolds  analogy  factor  rises  with  increas¬ 
ing  free-stream  turbulence  level.  It  should  be  pointed  out,  however,  that  the 
results  of  Refs.  2  and  14  clearly  show  that  the  effects  of  free-stream  turbulence 
on  turbulent  boundary  layers  are  not  dependent  on  turbulence  intensity  alone.  For  a 
fixed  free-stream  intensity  the  largest  impact  on  a  turbulent  boundary  layer  results 
if  the  integral  scale  of  the  turbulence  is  about  equal  to  the  boundary  layer  thick¬ 
ness.  Turbulence  with  integral  scales  significantly  smaller  or  larger  than  the 
boundary  layer  thickness  will  produce  reduced  effects. 
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The  experimental  data  discussed  previously  in  this  report  was  used  to  assess  the 
capability  of  a  boundary  layer  analysis  for  predicting  the  effect  of  free-stream  tur¬ 
bulence  on  momentum  and  thermal  boundary  layers.  Previously,  Blair  and  Werle  (Ref.  2) 
examined  the  effects  of  free-stream  turbulence  on  zero  pressure  gradient  flows.  They 
also  evaluated  the  ability  of  a  finite  difference  code  (Ref.  28),  which  used  a  tur¬ 
bulence  model  of  McDonald  et  al .  (Refs.  29  and  30),  to  predict  surface  heating  and 
skin  friction.  The  present  analytical  investigation,  which  is  a  continuation  of  the 
work  initiated  by  Blair  and  Werle,  makes  use  of  a  boundary  layer  analysis  (ABLE  - 
Analysis  of  the  Boundary  Layer  Equations)  recently  developed  by  Edwards,  Carter  and 
Werle  (Ref.  31).  This  new  boundary  layer  analysis  contains  the  McDonald  et  al.  tur¬ 
bulence  model  (Refs.  29  and  30)  utilized  by  Blair  and  Werle's  (Ref.  2)  previous  work. 

In  addition,  it  was  demonstrated  in  Ref.  31  that  results  obtained  from  the  ABLE 
analysis  and  the  boundry  layer  procedure  employed  by  Blair  and  Werle  were  in  excellent 
agreement  for  zero  pressure  gradient  flows.  In  the  present  study,  the  capability  of 
the  ABLE  code  to  accurately  predict  mean  flow  velocity,  mean  flow  temperature,  Reynolds 
shear  stress  and  turbulent  heat  transport  profiles  is  determined.  In  addition,  the 
turbulent  Prandtl  number  distribution  deduced  from  the  experimental  measurements  dis¬ 
cussed  earlier  in  this  report  is  used  in  the  boundary  layer  analysis  and  its  effect  on 
surface  heating  is  evaluated. 


1.  Prediction  Method 

The  ABLE  boundary-layer  code  provides  a  rapid  computation  of  two  dimensional  or 
axi symmetric  boundary-layer  flows  subject  to  a  prescribed  distribution  of  edge  Mach 
number,  streamwise  velocity,  or  static  pressure.  At  the  surface  a  distribtuion  of 
either  wall  temperature  or  heat  flux  may  be  imposed.  This  analysis  is  applicable  to 
attached  flows  which  are  laminar,  transitional,  or  turbulent.  A  detailed  descrip¬ 
tion  of  the  theory  used  in  the  ABLE  code  is  given  in  Ref.  31  and  a  flow  chart  of  the 
code  is  shown  in  Fig.  13.  An  implicit  finite-difference  technique  is  used  in  the 
ABLE  code  to  solve  the  boundary-layer  equations  which  are  written  in  nondimensional 
form  for  two  dimensional  flow  as  follows. 
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In  the  above  equations,  s  is  the  coordinate  along  the  surface,  n  is  the  coordinate 
normal  to  the  surface,  u  is  the  streamwise  velocity,  v  is  the  normal  velocity,  o  is 
the  static  density,  P  is  the  static  pressure,  h  is  the  static  enthalpy  and  H  is  the 
total  enthalpy  where  (in  non-dimensional  form) 


H:  t  +  j  u2  h  +  j  u2 


The  placement  of  a  bar  over  several  terms  is  used  to  denote  the  time  average  of 
various  turbulent  fluctuating  quantities  which  are  generally  considered  to  represent 
the  dominant  Reynolds  stress  terms  in  the  turbulent  boundary  layer  equations. 

The  ABLE  code  currently  contains  two  turbulence  models,  the  Cebeci-Smith  alge¬ 
braic  model  (Ref.  32)  and  the  McDonald  et  al.,  one  equation  turbulence  model  (Refs. 

29  and  30).  Both  models  are  based  on  an  eddy  viscosity  concept  in  which  the  Reynolds 
shear  stress  is  related  to  the  mean  flow  velocity  gradient  by 
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In  addition,  the  turbulent  heat  transport  is  related  to  the  Reynolds  shear 
stress  and  mean  flow  quantities  through  a  Reynolds  analogy  type  of  argument 
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where  Prt  is  the  turbulent  Prandtl  number.  This  code  presently  contains  two  transi¬ 
tion  models,  the  first  of  which  is  the  Dhawan  and  Narasimha  (Ref.  33)  forced  transi¬ 
tion  model  which  requires  the  specification  of  the  start  and  length  of  transition. 
The  second  model  is  a  natural  transition  model  developed  by  McDonald  and  Fish  (Ref.. 
29)  where  the  prediction  of  transition  is  controlled  by  the  integrated  form  of  the 
turbulence  kinetic  energy  equation.  In  the  present  investigation  the  ABLE  code  is 
applied  to  the  experimental  flows  discussed  previously  using  the  one  equation  turbu¬ 
lence  model  of  McDonald  and  Kreskovsky  (Ref.  30),  the  details  of  which  are  given  in 
the  next  section. 


2.  Turbulence  Model 

The  one  equation  turbulence  model  of  McDonald  and  Kreskovsky  (Ref.  30)  permits 
the  effect  of  free-stream  turbulence  to  be  included  in  the  computed  boundary  layer 
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analysis.  This  turbulence  model,  which  accounts  for  the  proper  approach  of  the  tur¬ 
bulence  le -e-'  in  the  outer  region  of  the  boundary  layer  to  the  local  edge  value,  is 
based  on  the  integral  form  of  the  turbulence  kinetic  energy  equation.  This  model  is 
an  extension  of  a  turbulence  model  developed  earlier  by  McDonald  and  Fish  (Ref.  29). 
The  eddy  viscosity  coefficient  is  expressed  in  nondimensional  form  as 

<£>=  f  ' >*£  (23) 

where  Re  is  a  reference  Reynolds  number,  and  £,  the  local  mixing  length,  is  expressed 
as  a  function  of  the  mixing  length  at  the  boundary  layer  edge  through  the  relation 

{tQnh(l7)  (‘-tonh  (l^))['-cosf8r)]}  (24) 
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where  <  is  the  von  Karman  constant,  6  is  the  boundary-layer  thickness,  and  is  the 
"shear  stress"  thickness  which  is  defined  as  the  first  location  from  the  outer  edge 
of  the  boundary  layer  where 
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and  Tmax  is  the  maximum  shear  stress  at  each  streamwise  location. 

The  local  value  of  is  obtained  through  the  solution  of  the  integral  form  of 
the  turbulence  kinetic  energy  equation  which  is  expressed  in  nondimensional  form  as 
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where  the  subscript  e  denotes  the  flow  quantities  at  the  boundary-layer  edge.  The 
influence  of  the  free-stream  turbulence  in  the  one  equation  model  comes  through  the 
term,  qe?,  which  acts  as  a  source  term  in  Eq.  27.  In  the  above  relation,  L  is  a 
dissipation  length  and  a],  a2,  and  83  are  structural  coefficients  that  relate  the 
Reynolds  shear  stress  and  turbulence  intensity  components  to  the  turbulent  kinetic 
energy  as  suggested  by  Townsend  (Ref.  34)  and  Bradshaw  and  Ferris  (Ref.  19).  These 
coefficients  are  given  by 
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where  McDonald  and  Fish  (Ref.  29)  suggest  the  values  a2 
be  expressed  as  'cl lows 

°o  (w) 

°i  - - 7TI - : 

,  +  0o(foo“‘) 

where 

a0  =  0.0115 

and  Re,  which  is  referred  to  as  the  Reynolds  number  based  on  momentum  thickness  by 
Shamroth  and  McDonald  (Ref.  35)  but  is  in  fact  a  correlation  given  in  terms  of  an 
integrated  turbulent  Reynolds  number  R-p. 
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where  1S  an  estimate  of  the  inner  wall  layer.  The  computational  transition  pro¬ 
cess  is  controlled  by  the  structural  coefficient  ai  as  it  varies  from  zero  in  laminar 
flow  to  .15  in  in  fully  turbulent  flow.  For  the  present  analysis,  the  structural 
coefficients  a],  and  a3  are  constant  over  the  boundary  layer  using  the  values 
suggested  by  McDonald  and  Fish  (Ref.  29).  However  the  experimental  structural  coeffi¬ 
cients  were  observed  to  vary  across  the  boundary  layer  (Figs.  8(a)  and  8(b)  and  the 
effect  of  varying  the  structural  coefficients  in  the  analysis  should  be  assessed  in 
the  future.  A  detailed  description  of  the  one  equation  turbulence  model  is  given  in 
Ref.  36. 


3.  Turbulent  Prandtl  Number  Model 

In  the  present  investigation  an  evaluation  is  made  of  the  effect  of  a  variable 
turbulent  Prandtl  number,  Prt,  across  the  boundary  layer  on  the  ABLE  cod?  prediction 
of  the  turbulent  heat  transport  and  surface  heating.  Three  different  functional 
forms  of  Prt  have  been  applied  in  this  investigation.  They  are  given  by 

1)  McDonald  (Ref.  28) 

,P(2ip),P(-!^p)  (42, 
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where 


P(x>=  -^[i  *  erf(-5=)] 


2)  Rotta  (Ref.  22) 


Pf\ -  .95  -  45(n/S)2 


3)  Present  experimental  investigation 


(43) 


2 

Prt-[  95  -  45(n/8)2  ][i  +  2Tu]  -  5Tu/cosh[ion/s] 


(44) 


where 


Tu=  i/3qe2 


A  comparison  of  these  three  turbulent  Prandtl  number  distributions  are  shown  in 
Fig.  14.  McDonald's  function  (Eq.  42)  has  a  maximum  value  of  1.7  at  the  wall  and 
decreases  rapidly  in  the  laminar  sublayer  of  the  turbulent  boundary  layer  but  is 
nearly  constant  and  equal  to  .9  for  n+  >  50.  Rotta's  function  (Eq.  43)  has  a  maximum 
value  of  Prt  of  .95  at  the  wall  and  decreases  linearly  with  respect  of  (n/5)2  to 
the  edge  of  the  boundary  layer  to  a  value  of  .5.  The  Pr^  distribution  obtained  in 
the  experimental  portion  of  the  present  investigation  (Eq.  44)  is  essentially  a  modi¬ 
fication  of  Rotta's  distribution  to  account  for  the  effects  of  free-stream  turbulence. 
An  assessment  of  the  accuracy  of  each  of  these  turbulent  Prandtl  number  formulations 
is  made  in  the  next  section  using  the  ABLE  code  and  the  experimental  results  presented 
above.  However,  since  the  turbulent  heat  transport  h'v'  is  modeled  in  terms  of  u'v', 
the  turbulent  Prandtl  number  and  the  normal  derivatives  of  u  and  h  as  given  in  Eq. 

22,  then  all  of  these  quantities  are  compared  with  the  experimental  data  before  an 
assessment  is  made  of  the  effect  of  Pr^  on  the  calculation  of  the  turbulent  heat 
transport. 


DISCUSSION  OF  ANALYTICAL  RESULTS 


A  series  of  calculations  have  been  made  with  the  ABLE  code  for  the  flow  over  a 
heated  flat  plate  for  each  of  the  nominal  inlet  free-stream  turbulence  levels  of  1%, 
2%,  4%  and  6 %  generated  by  the  use  of  inlot  turbulence  Grids  1,  2,  3  and  4,  respec¬ 
tively.  A  calculation  was  not  performed  for  the  case  of  the  flow  with  .25  percent 
inlet  turbulence  since  it  was  concluded  by  Blair  and  Werle  (Ref.  2)  that  the  transi¬ 
tion  process  of  this  flow  is  three  dimensional  and  thus  the  turbulence  model  of 
McDonald  et  al.  (Refs.  29  and  30)  cannot  accurately  predict  the  location  and  length 
of  the  transition  region.  For  each  case,  a  calculation  is  made  with  the  ABlE  code 
for  each  of  the  turbulent  Prandtl  number  formulations  discussed  in  the  previous 
section.  The  mean  flow  quantities,  Reynolds  shear  stress  and  turbulence  kinetic 
energy  predicted  from  the  ABLE  code  were  found  to  be  insensitive  to  the  different 
turbulent  Prandtl  number  formulations;  hence,  these  quantities  are  presently  for  only 
the  present  Pr^  formulation  given  in  Eq.  44.  This  result  was  expected  since  the 
experimental  flows  are  low  speed  and  thus  the  momentum  equation  (Eq.  19)  is  essen¬ 
tially  uncoupled  from  the  energy  equation  (Eq.  20). 

For  all  of  the  test  cases  analyzed  in  this  investigation  with  the  ABLE  code,  the 
gas  is  assumed  to  be  air  with  a  constant  ratio  of  specific  heats,  v,  equal  to  1.4  and 
a  constant  Prandtl  number  equal  to  .72.  The  von  (Carman  constant  for  turbulent  flow 
is  set  to  .43  as  suggested  by  McDonald  and  Kreskovsky  (Ref.  30).  The  following  flow 
conditions  were  used  in  all  test  cases 


Ue  *  100  ft/ sec 
PTe  -  14.78  lb/in2 


and  the  streamwise  variation  of  the  free-stream  turbulence  for  each  flow  with  a  spe 
cified  inlet  grid  is  obtained  from  the  following  expression 


Tu=  78  ( 


2.54x  +132  ,~5/7 
b  ' 


(45) 


where  b  (grid  bar  width)  =  .48,  1.27,  3.81,  and  5.08  for  Grids  1,  2,  3,  and  4,  respec¬ 
tively.  This  relation  was  shown  to  be  accurate  in  earlier  testing  reported  in  Ref.  2. 
The  measured  wall  temperature  levels  are  tabulated  in  Table  I  of  Ref.  2,  with  the 
free-stream  static  temperature  set  to  530°R  (TTe  =  530.83°R)  for  all  calculations. 

The  temperature  distributions  were  numerically  smoothed  to  eliminate  spurious  varia¬ 
tions  in  the  computed  wall  results  due  to  minor  experimental  error.  The  smoothed 
temperature  distributions  were  used  as  input  to  the  ABLE  code.  Comparison  of  the 
measured  and  smoothed  temperature  distributions  for  each  of  the  flow  cases  are  shown 
in  Fig.  15.  The  smoothing  procedure  is  a  least  squares  polynomial  curve  fit  described 
in  Ref.  37. 

A  computational  mesh  consisting  of  101  grid  points  in  the  normal  direction  and 
100  points  in  the  streamwise  direction  was  used  in  each  of  the  calcualtions.  A  grid 
stretching  based  on  a  geometric  progression  was  applied  in  each  direction  to  insure 
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that  a  fine  grid  distribution  was  placed  in  the  high  gradient  regions.  The  initial 
profile  for  the  boundary  layer  calculation  is  the  Blasius  profile  (Ref.  38)  which 
was  imposed  at  the  flat  plate  leading  edge.  In  addition,  two  iterations  per  stream- 
wise  station  are  applied  in  the  computational  procedure  due  to  large  streamwise 
temperature  gradients  which  are  encountered  in  the  transition  region  of  the  flow.  The 
current  calculations  required  approximately  1  minute  of  CPU  time  on  a  UNIVAC  1180 
operating  system  to  compute  the  flow  over  the  8  foot  length  of  the  test  section. 

In  this  section  comparisons  are  presented  in  Figs.  16-27  between  the  results 
obtained  from  the  ABLE  code  with  those  measured  experimentally  both  in  the  present 
investigation  and  in  the  previous  investigation  by  Blair  and  Werle  (Ref.  2).  These 
comparisons  are  made  for  zero  pressure  gradient  flows  for  the  two  cases  with  inlet 
turbulence  levels  of  2  percent  and  6  percent.  The  following  quantities  are  compared: 

1.  skin  friction 

2.  displacement  thickness 

3.  momentum  thickness 

4.  mean  velocity  profile 

5.  mean  temperature  profile 

6.  Reynolds  shear  stress  profile 

7.  turbulent  kinetic  energy  profile 

8.  profiles  of  the  components  of  turbulence  intensity 

9.  turbulent  Prandtl  number  profile 

10.  turbulent  heat  transport  profile 

The  profile  comparisons  are  shown  only  at  X  =  68  inches  from  the  leading  edge  of  the 
test  section;  similar  comparisons  were  obtained  at  the  other  measuring  stations. 

Figure  16  is  a  comparison  of  the  computed  skin  friction  coefficient  distribution 
with  that  obtained  experimentally  along  the  flat  plate  surface  for  inlet  turbulence 
levels  of  2  percent  and  6  percent.  The  theoretical  distribution  is  slightly  higher 
(approximately  5%)  than  the  experimental  distribution  for  the  flow  with  an  inlet  tur¬ 
bulence  of  2  percent  while  the  result  obtained  for  the  flow  with  an  inlet  turbulence 
of  6  percent  is  in  excellent  agreement  with  the  experimental  result.  A  comparison 
of  the  theoretical  and  experimental  displacement  and  momentum  thickness  distributions 
are  shown  in  Figs.  17  and  18,  respectively.  It  is  apparent  from  these  figures  that 
the  computed  integral  quantities  are  in  good  agreement  with  the  experiment  over  most 
of  the  test  section.  Figures  19  and  20  are  comparisons  of  the  computed  mean  velo¬ 
city  and  temperature  profiles  with  the  experimental  measurements  for  inlet  turbu¬ 
lence  levels  of  2  percent  and  6  percent.  In  all  profile  comparisons,  the  theoretical 
boundary  layer  thickness,  6,  was  used  to  nondimensionalize  the  normal  coordinate. 

It  is  observed  from  Figs.  19  and  20  that  the  computed  results  are  in  excellent 
agreement  with  the  experimental  results  for  both  flows.  From  Figs.  16-20  it  is  con¬ 
cluded  that  the  ABLE  boundary  layer  analysis  with  the  McDonald  and  Kreskovsky  turbu¬ 
lence  model  (Ref.  30)  produces  good  agreement  with  the  mean  flow  quantities  of  zero 
pressure  gradient  flows  with  various  levels  of  inlet  turbulence. 

In  order  to  determine  how  well  the  boundary  layer  analysis  will  predict  turbu¬ 
lent  fluctuating  quantities,  comparisons  of  computed  profiles  with  experimental  data 
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were  made  of  the  flows  wiJ  »  inlet  turbulence  levels  of  2  percent  and  6  percent. 

Figure  21  is  a  comparison  of  the  computed  Reynolds  shear  stress  profile  with  that 
obtained  experimentally.  Several  interesting  features  are  noted  in  this  figure.  By 
using  the  turbulence  model ing(Eqs.  28-32)  developed  by  McDonald  and  Kreskovsky  in  the 
computational  procedure,  the  predicted  Reynolds  shear  stress  in  the  inner  layer  of 
the  turbulent  boundary  layer  is  in  excellent  agreement  with  the  experimental  measure¬ 
ments.  The  theoretical  results  exhibit  the  same  shape  and  level  as  the  experimental 
results.  However,  in  the  outer  layer,  as  the  free-stream  turbulence  increases,  the 
deviation  between  the  computed  and  experimental  results  grows  significantly  as  it  is 
observed  that  the  computed  u'v'  0  whereas  the  experimental  u'v'  approaches  a  finite 
value.  Figure  22  is  a  comparison  of  the  measured  and  computed  turbulence  kinetic 
energy  distribtuions.  It  is  observed  that  the  computed  results  are  in  good  agreement 
with  the  experimental  data  for  both  levels  of  free-stream  turbulence.  This  result 
indicates  that  McDonald  and  Kreskovsky' s  modeling  of  free-stream  turbulence  in  the 
turbulence  kinetic  energy  equation  (Eq.  27)  captures  the  correct  shape  of  the  turbu¬ 
lence  kinetic  energy  across  the  boundary  layer  for  zero  pressure  gradient  flows  with 
different  levels  of  free-stream  turbulence.  Figures  23-25  are  comparisons  of  the 
measured  and  computed  components  of  turbulence  intensity,  u'u1,  v'v',  w‘w’.  From 
these  figures  it  is  observed  that  the  predicted  results  are  in  fair  agreement  with 
the  experimental  data  across  the  boundary  layer. 

An  anomaly  appears  in  the  results  shown  in  Figs.  19-25.  First,  from  Figs.  22-25, 
the  computed  turbulence  kinetic  energy  and  its  various  components  are  in  realtively 
good  agreement  with  the  experimental  results.  Secondly,  the  computed  mean  velocity 
(Fig.  19)  is  in  excellent  agreement  with  the  experiment.  However,  the  computed 
Reynolds  shear  stress  (Figs.  21)  has  a  significant  deviation  from  the  experimental 
results  in  the  outer  region  of  the  boundary  layer  as  the  inlet  turbulence  level 
increases.  This  difference  tends  to  suggest  that  the  eddy  viscosity  modeling  of  tur¬ 
bulence  (Eq.  21)  does  not  properly  model  flows  with  significant  levels  of  free-stream 
turbulence  since  by  applying  this  model,  the  computed  Reynolds  shear  stress  is  forced 
to  zero  at  the  edge  of  the  boundary  layer.  The  inability  of  the  turbulence  model 
to  predict  accurate  Reynolds  shear  stress  distributions  over  the  entire  boundary 
layer  will  affect  the  transition  model  (Eq.  39)  since  the  transition  process  is  con¬ 
trolled  by  an  integrated  form  of  the  turbulent  Reynolds  number  Rj  (Eq.  41).  An  in¬ 
vestigation  is  needed  to  determine  an  analytical  turbulence  model  that  will  properly 
represent  the  Reynolds  shear  stress  in  the  outer  region  of  the  turbulent  boundary 
layer  for  flows  with  significant  levels  of  free-stream  turbulence. 

A  comparison  of  McDonald's,  Rotta's,  and  the  present  turbulent  Prandtl  number 
formulation  with  the  experimentally  determined  distribution  is  shown  in  Fig.  26(a) 
and  Fig.  26(b)  for  flows  with  2  percent  and  6  percent  turbulence  levels,  respectively. 
In  Fig.  26(a)  it  is  observed  that  McDonald's  distribution,  which  is  essentially 
constant  over  the  boundary-layer  overpredicts  the  experimental  results  over  most  of 
the  boundary  layer  while  Rotta's  and  the  present  distribtuion  exhibit  the  same 
general  shape  and  level  as  the  experiment  distributions.  For  the  high  inlet  turbu¬ 
lent  flow  (Fig.  26(b)),  the  present  Pr^  distribution  exhibits  the  same  shape  and 
level  as  the  experimental  distribution  while  the  relatively  constant  McDonald 
distribtuion  does  not  have  the  same  shape  or  level  as  the  experimental  results. 
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Also,  Rotta's  distribution  exhibits  the  best  shape  of  the  experimental  data  but  is 
not  at  the  same  level.  Figures  27(a)  and  26(b)  are  comparisons  of  the  measured  and 
computed  turbulent  heat  flux  distributions  for  the  flows  with  turublence  levels  of 
2  percent  and  6  percent,  respectively.  From  these  figures,  it  is  noted  that  except 
near  the  wall,  the  computed  turbulent  heat  flux  distribution  across  the  boundary 
layer  is  essentially  the  same  for  each  of  the  three  different  Pr^  formulations  which 
are  used.  Figures  27(a)  and  27(b)  show  that  the  computed  results  are  in  reasonable 
agreement  with  the  experimental  results  except  at  the  edge  of  the  boundary  layer 
where  the  theoretical  results  go  to  zero  and  the  experimental  data  does  not.  Since 
the  turbulent  heat  flux  is  determined  from  Eq.  (22)  and  all  the  computed  quantities 
in  that  equation  are  in  reasonable  agreement  with  experimental  data  except  the 
Reynolds  shear  stress,  this  suggests  that  the  inaccuracy  of  predicting  the  turbuent 
heat  flux  in  the  outer  portion  of  the  boundary  layer  is  due  to  the  modeling  of  the 
Reynolds  shear  stress  in  the  outer  region  of  a  turbulent  boundary  layer  with  signi¬ 
ficant  levels  of  free-stream  turbulence.  Further  investigation  of  this  feature  of 
the  flow  is  needed. 

In  Figure  28,  comparisons  are  presented  between  measured  and  computed  Stanton 
number  distributions  for  the  flows  with  inlet  turbulence  levels  of  1,  2,  4  and  6 
percent.  The  following  observations  about  the  Stanton  number  prediction  in  the  fully 
turbulent  region  of  the  flows  are  made.  The  predicted  Stanton  number  using  McDonald's 
formulation  is  in  good  agreement  with  the  experiment  for  the  flow  with  an  inlet  tur¬ 
bulence  level  of  1  percent;  however,  as  the  inlet  turbulence  level  increases,  the 
computational  procedure  using  McDonald's  formulation  underpredicts  the  measured 
Stanton  number  distribution.  The  computational  procedure  using  Rotta's  formulation 
overpredicts  the  Stanton  number  for  the  flow  with  the  1  percent  inlet  turbulence 
level.  However,  as  the  inlet  turbulence  level  increases  the  computed  results  tend 
to  slightly  overpredict  the  measured  Stanton  number  distribution.  The  present 
turbulent  Prandtl  number  formulation  yields  essentially  the  same  results  as  the 
computed  with  Rotta's  formulation  since  the  present  formulation  is  a  perturbation 
of  Rotta's  Prt  formulation.  The  results  shown  in  Fig.  28  indicate  that  the  computa¬ 
tion  using  the  present  distribution  shows  no  marked  improvement  over  the  computation 
using  Rotta's  distribution.  However,  for  the  flows  with  inlet  turbulence  of  2  per¬ 
cent  or  larger,  computations  using  either  Rotta's  or  the  present  formulation  result 
in  predicted  Stanton  number  distributions  which  are  in  better  agreement  with  the 
experimental  results  than  that  obtained  using  McDonald's  formulation.  A  further 
indication  of  the  advantage  of  Rotta's  or  the  present  formulation  is  shown  in  Fig.  29 
where  the  predicted  Reynolds  analogy  factor,  2  St/Cf,  from  calculations  using  the 
present  and  McDonald's  turbulent  Prandtl  number  distributions  are  compared  to  the 
experimentally  deduced  Reynolds  analogy.  In  this  figure  it  is  observed  that  the  cal¬ 
culation  using  the  present  (or  Rotta's)  formulation  predicted  Reynolds  analogy 
factors  that  are  in  better  agreement  with  the  experiment  than  the  computation  which 
uses  McDonald's  formulation.  The  overall  implication  of  these  results  is  that  the 
analysis  using  the  McDonald  and  Kreskovsky  turbulence  model  (Ref.  30)  with  either 
Rotta's  or  the  present  turbulent  Prandtl  number  formulation  can  accurately  repre¬ 
sent  the  momentum  and  energy  transport  mechanisms  for  zero  pressure  gradient  flows 
in  the  wall  region  of  the  boundary  layer  but  that  there  is  a  severe  weakness  in  its 
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ability  to  represent  the  momen'1  j.i  and  energy  transport  mechanisms  in  the  region  near 
the  edge  of  the  boundary  layer  .or  flows  with  significant  levels  of  free-stream 
turbulence. 
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The  present  program  was  designed  to  examine,  both  experimentally  and  analyti¬ 
cally,  the  effect  of  free-stream  turbulence  on  the  heat  transfer  through  turbulent 
boundary  layers.  The  experimental  test  conditions  for  the  present  program  were 
intended  to  reproduce  cases  for  which  numerous  other  experimental  data  had  been  ob¬ 
tained  under  an  earlier  AFOSR  contract  (Ref.  2).  Measurements  of  multi -component 
free-stream  turbulence  intensities,  test  surface  Stanton  number  distributions,  tran¬ 
sition  Reynolds  numbers  and  boundary  layer  integral  thicknesses  were  in  excellent 
agreement  with  the  respective  quantities  of  the  earlier  contract.  It  has  been  con¬ 
cluded  that  these  present  measurements  can  be  viewed  as  additional  data  for  the  same 
test  conditions  as  were  previously  studied.  A  number  of  comparisons  were  made  be¬ 
tween  low  free-stream  turbulence  boundary  layer  turbulence  data  obtained  in  the 
present  study  and  similar  results  from  other  investigations.  These  comparisons 
showed  excellent  agreement  indicating  that  the  present  boundary  layer  turbulence 
data  are  of  high  quality. 

The  conclusions  reached  from  the  experimental  measurements  obtained  for  higher 
levels  of  free-stream  turbulence  were  as  follows: 

1.  The  present  data  indicate  a  progressive  increase  of  boundary  layer 
turbulence  kinetic  energy  with  increasing  free-stream  turbulence. 

Increased  levels  of  turbulence  kinetic  energy  were  measured  across 
the  entire  thickness  of  the  boundary  layer.  These  results  are  in 
agreement  with  data  from  other  independent  studies. 

2.  Both  the  u'  and  w '  components  of  turbulence  increased  progressively 
with  increasing  free-stream  turbulence  level.  The  u1  component 
increased  more  than  the  w'  component.  The  vertical  component  (v1), 
however,  was  essentially  constant  and  independent  of  free-stream 
turbulence  level  for  the  inner  half  of  the  boundary  layer. 

3.  The  ratio  of  shear  stress  to  turbulence  kinetic  energy  decreased 
across  the  entire  boundary  layer  with  increasing  free-stream  turbu¬ 
lence  level.  The  decrease  was  most  extreme  over  the  outer  60  percent 
of  the  boundary  layer. 

4.  The  effects  of  free-stream  turbulence  level  on  the  ratios  of  the 
direct  stress  components  to  the  turbulence  kinetic  energy  were  to 
a)  increase  u'u’/q^,  b)  decrease  v'v'/ q2  and  c)  leave  wV/qZ~ 
nearly  constant. 

5.  Reynolds  stress  distribution  measurements  indicated  that  at  high 
levels  of  free-stream  turbulence  the  turbulent  shear  stresses  extend 
beyond  the  mean  velocity  boundary  layer.  The  present  data  and 
results  from  other  sources  indicate  an  increase  in  turbulent  shear 
at  the  boundary  layer  edge  directly  proportional  to  the  free-stream 
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turbulence  level.  Flatness  factor  measurements  indicated  that  as 
the  free-stream  turbulence  level  was  increased  the  "border"  between 
the  fuid  in  the  boundary  layer  and  the  free-stream  fluid  became 
less  distinct. 

6.  Measurements  of  the  boundary  layer  turbulent  Prandtl  number  distri¬ 
bution  for  the  case  of  the  low  free-stream  turbulence  were  in  good 
agreement  with  a  model  suggested  by  Rotta.  The  present  data  indicate 
that  as  the  free-stream  turbulence  level  was  increased,  the  near-wall 
Prt  decreased  while  Prt  over  the  outer  region  of  the  boundary  layer 
slightly  increased.  A  correlation,  Prt  (y/6,  Te),  which  fit  the  ob¬ 
served  data  reasonably  well  was  suggested. 

The  experimental  data  was  used  to  assess  the  capability  of  a  boundary-layer  com¬ 
puter  program,  ABLE  (Analysis  of  the  Boundary  L.ayer  Equations)  for  predicting  the 
effect  of  free-stream  turbulence  on  momentum  and  thermal  boundary  layers.  In  addi¬ 
tion  the  turbulent  Prandtl  number  formulation  deduced  from  the  experimental  measure¬ 
ments  was  used  in  the  boundary  layer  analysis  and  its  effect  on  surface  heating  was 
determined.  The  following  conclusions  were  reached  from  the  theoretical  portion  of 
this  investigation: 

1.  The  modeling  of  free-stream  turbulence  in  the  one  equation  turbu¬ 
lence  model  of  McDonald  and  Kreskovsky  captures  the  correct  shape 
and  level  of  the  turbulence  kinetic  energy. 

2.  For  increased  levels  of  free-stream  turbulence,  the  Reynolds  shear 
stress  and  turbulent  heat  flux  determined  from  the  turbulence  model 
is  signif icantly  smaller  than  that  observed  experimentally  in  the 
wake  region  of  the  turbulent  boundary  layer.  This  discrepancy 
could  be  due  to  the  eddy  viscosity  concept  used  in  McDonald  and 
Kreskovsy's  model  and  further  investigation  of  turbulence  models  is 
needed. 

3.  Analytical  calculations  using  either  Rotta 's  turbulent  Prandtl 
number  correlation  or  the  correlation  of  the  present  investigation 
predicted  Reynolds  analogy  factors  (2St/Cf)  that  are  in  reasonable 
agreement  with  experimental  measurements  and  accurately  predict  the 
increase  in  surface  heat  transfer  due  to  increased  free-stream 
turbulence. 
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LIST  OF  SYMBOLS 


a-| ,  &2t  83 ,  84,  a-j  ^ a 

CP 

d 

Fuv 

Fuh 

Fvv 

Fwh 

Fj 

G1  ,e 
h 

h 

h '  v 1 
H 

l 

£ 

L 

n 

n+ 

Nu 

P 

Prt 

q2 

r 

Re 

Ree 

«T 

Re 


Turbulence  structural  coefficients 

Specific  heat  at  constant  pressure 

Hot  wire  sensor  diameter 

Wall  damping  function  in  turbulence 

Flatness  factor  (vertical  x  probe)  =  u'V(u'2)^ 

Flatness  factor  (horizontal  x  probe)  =  u'4/(u'2)2 

Flatness  factor  (vertical  x  probe)  =  v'4/(v'2)2 

Flatness  factor  (horizontal  x  probe)  =  w'4/(w'2)2 

Flatness  factor  (tri-x  probe)  =  t"fy(t'2)2 

Turbulence  structural  coefficient  (Eq.  11) 

Height  of  hot  wire  sensor  array 
Static  enthalpy 
Reynolds  thermal  flux 
Total  enthalpy 

Active  length  of  hot  wire  sensor 
Mixing  length  function  in  turbulence  model 
Free-stream  mixing  length 
Dissipation  length  scale 
Distance  normal  to  surface 
Dimensionless  normal  distance  to  surface, 

Nusselt  number  of  hot  wire  sensor 
Static  pressure 

Turbulent  Prandtl  number  (Eq.  9) 

Heat  flux 

Turbulence  kinetic  energy 

Separation  distance  between  hot  wire  sensors 

Reynolds  number 

Reynolds  number  based  on  momentum  thickness 
Turbulent  Reynolds  number 
McDonald's  correlation  of  Ree 
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u * u * ,  v'v1 ,  w'w‘ 


u ' ,  v' ,  w' 


-u’v’ 


Coordinate  along  v\ll  surface 
Skewness  factor  (vertical  x  probe)  = 

Skewness  factor  (horizontal  x  probe)  =  u '3/(u ,2)^y/^ 
Skewness  factor  (vertical  x  probe)  =  v‘3/(v'2)^//^ 

<  ■  ■  «3  /rt 

Skewness  factor  (horizontal  x  probe)  =  w'3/(w'2)J/ 
Skewness  factor  (tri-x  probe)  =  t*3/(t'2)^ 

Mean  static  temperature 
Free-stream  turbulence  level 
Total  temperature 
Friction  temperature  =  qw/pwCpUT 
Dimensionless  temperature  =  Tw-T/Tt 
Fluctuating  temperature 
Streamwise  velocity 
Mean  streamwise  velocity 
Friction  velocity 
Dimensionless  velocity  -  U/UT 
Components  of  turbulent  intensity 

Streamwise  normal  and  transverse  fluctuating  velocities 
Reynolds  shear  stress 
Normal  velocity 

Transport  velocity  of  turbulence  kinetic  energy  (Eq.  8) 
Transport  velocity  of  turbulent  shear  stress  (Eq.  7) 
Distance  from  plate  leading  edge 
Distance  from  wall 

Dimensionless  distance  from  wall  =  yllT/v 

Boundary  layer  thickness 

Thermal  boundary  layer  thickness 

Thickness  of  inner  wall  region  of  boundary  layer 

Boundary  layer  displacement  thickness 

Boundary  layer  momentum  thickness 

von  Karman  constant  for  velocity  of  law-of-the-wall 

von  Karman  constant  for  temperature  law-of-the-wall 


Subscripts 


Molecular  viscosity 
Eddy  viscosity 
Kinematic  viscosity 

Wake  strength  for  velocity  boundary  layer 
Wake  strength  for  temperature  boundary  layer 
Density 
Shear  stress 

Correlation  coefficient  =  u'^v'/u'^  ■J& 
Correlation  coefficient  =  u'^w'/u'^ 
Correlation  coefficient  =  u 'v ' 2/u '2  fy'2 
Correlation  coefficient  =  u 'w'2/u '2  \/w 


Freestream 


0.995 


where  U  =  0.995  Ue 
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HEAT  EXCHANGER 


Figure  1.  United  Technologies  Research  Center  Boundary  Layer  Wind  Tunnel 


Figure  2.  Photograph  of  UTRC  Boundary  Layer 
Wind  Tunnel  Test  Section 


RIGID  URETHANE  FOAM 


Figure  3.  Instrumentation  Diagram  for  the  Uniform  Heat  Flux  Flat  Wall  Model 
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Figure  7.  Distribution  of  the  Components  of  Boundary  Layer  Turbulence 
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Figure  8A.  Distribution  of  the  Turbulence 
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Figure  13.  Flow  Chart  for  ABLE  Code  -  Module  MAIN 
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Figure  14.  Comparison  of  Turbulent  Prandtl  Number  Distributions  for  Flow 
with  2 %  Freestream  Turbulence 
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Figure  15.  Numerically  Smoothed  Experimental  Wall  Temperature  Distribution 
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a)  FLOW  WITH  6%  INLET  TURBULENCE 


b)  FLOW  WITH  2%  INLET  TURBULENCE 


Figure  16.  Comparison  of  Theoretical  Skin  Friction  with  Experimental  Data 
for  Different  Inlet  Turbulence  Levels 
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Figure  17.  Comparison  of  Theoretical  Displacement  Thickness  with  Experimental 
Data  for  Different  Inlet  Turbulence  Levels 
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Figure  18.  Comparison  of  Theoretical  Momentum  Thickness  with  Experimental 
Data  for  Different  Inlet  Turbulence  Levels 
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Figure  19.  Comparison  of  Theoretical  Mean  Velocity  With  Experimental  Data 
at  X  =  68  inches  for  Different  Inlet  Turbulence  Levels 
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Figure  20.  Comparison  of  Theoretical  Mean  Temperature  with  Experimental 
Data  at  X  =  68  inches  for  Different  Inlet  Turbulence  Levels 
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Comparison  of  Theoretical  Reynolds  Shear  Stress  with  Experimental 
Data  at  X  =  68  inches  for  Different  Inlet  Turbulence  Levels 
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Figure  22.  Comparison  of  Theoretical  Turbulent  Kinetic  Energy  with  Experimental 
Data  at  X  =  68  inches  for  Different  Inlet  Turbulence  Levels 
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Figure  23.  Comparison  of  Theoretical  Component  of  Turbulent  Intensity,  u’v', 
with  Experimental  Data  at  X  =  68  inches  for  Different  Inlet 
Turbulence  Levels 
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Figure  24.  Comparison  of  Theoretical  Component  of  Turbulent  Intensity,  vV, 
with  Experimental  Data  at  X  =  68  inches  for  Different  Inlet 
Turbulence  Levels 
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Figure  25.  Comparison  of  Theoretical  Component  of  Turbulent  Intensity,  wV, 
with  Experimental  Data  at  X  =  68  inches  for  Different  Turbulence 
Levels 
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Comparison  of  Theoretical  Turbulent  Prandtl  Number  with  Experimental 
Data  at  X  =  68  inches  for  2%  Inlet  Turbulence  Level 
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Figure  27a.  Comparison  of  Theoretical  Thermal  Heat  Flux  with  Experimental  Data 
at  X  *  68  inches  for  a  2%  Inlet  Turbulence  Level 
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Comparison  of  Theoretical  Thermal  Heat  Flux  with  Experimental  Data 
at  X  *  68  inches  for  a  6%  Inlet  Turbulence  Level 
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Comparison  of  Theoretical  Stanton  Number  with  Experimental  Data 
for  2%  Inlet  Turbulence  Level 
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Figure  28c.  Comparison  of  Theoretical  Stanton  Number  with  Experimental  Data 
for  4%  Inlet  Turbulence  Level 
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Figure  28d.  Comparison  of  Theoretical  Stanton  Number  with  Experimental  Data 
for  6%  Inlet  Turbulence  Level 
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Figure  29.  Theoretical  Prediction  of  Reynolds  Analogy  Factor  -  Varying 
Edge  Turbulence  Levels 
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APPENDIX  A  -  ERROR  ANALYSIS 


In  this  program  mean  velocity  (U)  profiles  were  determined  by  means  of  Pitot, 
single  wire,  x  wire  and  3  sensor  hot  wire  probes.  Mean  temperatures  (T)  were  measured 
by  both  3  sensor  hot  wire  and  thermocouple  probes .  Fluctuating  velocities  and 
Reynolds  stress  distributions  (u1,  v‘,  w u'v1)  were  determined  both  with  x  wire  and 
3  wire  techniques. 

Assessment  of  absolute  errors  for  the  analog  signals  measured  in  this  program 
would  be  a  relatively  straightforward  matter.  For  example,  the  possible  errors  in 
measured  pressures  from  the  pitot  probes  or  the  recorded,  digitized  voltages  from  the 
hot-wire  anemometer  can  be  computed  from  the  individual  expected  uncertainties.  Com¬ 
putation  of  the  absolute  errors  of  the  measured  physical  quantities  (e.g.,  u',  v',  t ' ) , 
however,  is  practically  impossible  because  the  true  accuracy  of  factors  such  as  Pitot 
probe  wall  proximity  corrections,  Pitot  probe  turbulence  corrections,  hot  wire  wall 
radiation  effects,  high  turbulence  sensor  cross-talk,  etc.  are  unknown.  For  this 
reason  the  uncertainties  for  the  various  quantities  measured  for  this  program  will  be 
assessed  by  (1)  comparing  the  measured  quantities  with  independently  determined  or 
computed  results  or  (2)  comparisons  of  like  quantities  measured  using  different  probes 
and  instrumentation  techniques. 

Mean  velocity  and  temperature  profile  data  obtained  with  different  measurement 
techniques  (four  techniques  for  velocity,  two  techniques  for  temperature)  are  presented 
in  Fig.  A-l .  These  profile  data  were  all  obtained  at  X  =  84  inches  at  three  levels 
of  free-stream  turbulence.  In-depth  descriptions  of  the  Pitot  and  thermocouple  probes 
and  data  systems  used  for  the  mean  velocity  and  temperature  profiles  are  provided  in 
Ref.  2.  These  probes  were  designed  and  constructed  specifically  for  these  types  of 
boundary  layer  flows  and  a  number  of  well  established  near-wall  correction  terms  (see 
Ref.  2)  were  applied  to  the  data.  In  addition,  a  number  of  comparison  checks  (see 
Refs.  1  and  2)  showed  that  these  mean  profile  data  were  very  accurate.  For  these 
reasons  the  mean  velocity  profiles  from  the  Pitot  probe  and  the  mean  temperature 
profiles  from  the  thermocouple  probe  were  selected  as  the  "true"  respective  profiles. 
For  the  data  of  Fig.  A-l  all  the  other  velocity  profile  measurements  were  compared  to 
the  Pitot  probe  data  while  the  triple  sensor  temperature  data  were  compared  to  the 
temperature  profiles  from  the  thermocouple  probe. 

An  examination  of  the  mean  velocity  profile  data  of  Fig.  A-l  indicates  that  85 
percent  of  all  the  measurements  fell  within  ±  3  percent  of  the  "true"  Pitot  profile. 
Only  one  set  of  data,  the  vertical  x  probe  results  for  Te  =  0.2  percent,  had  any  dis¬ 
crepancies  larger  than  5  percent.  Discrepancies  in  the  mean  temperature  profile 
measurements  were  slightly  larger  than  those  for  the  velocity  measurements  with  only 
70  percent  of  the  measurements  falling  within  ±  5  percent  of  the  thermocouple  probe 
profile.  The  3  sensor  probe  temperature  data  also  showed  a  clear  bias  to  read 
slightly  lower  temperatures  than  the  thermocouple  probe.  The  relatively  larger  errors 
for  the  temperature  profiles  as  compared  to  the  velocity  profiles  is  not  surprising. 
Because  of  the  relatively  large  size  of  the  3  sensor  probes  they  span  a  gradient  of 
both  velocity  and  temperature  in  the  boundary  layer  flow.  The  data  reduction  system 
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is  forced  to  assume  that  a  single  effective  velocity  and  flow  temperature  apply  over 
the  entire  probe. 


Plots  of  the  average  bias  (xi-x)  and  standard  error  V(xi-x)2  (xi  is  the 
measured  quantity  and  x  is  the  "true")  quantity  at  a  given  profile  location)  for  the 
various  profiles  are  given  in  Fig.  A-2.  These  results  are  plotted  as  a  function  of 
the  ratio  between  the  overall  boundary  layer  thickness  and  the  probe  sensor  height 
(h).  The  overall  profile  errors  are  plotted  in  Fig.  A-2,  were  largest  for  the  thin¬ 
nest  boundary  layers  (with  the  relatively  steepest  velocity  and  temperature  gradients). 
In  addition,  the  largest  local  bias  errors  of  Fig.  A-l  were  located  near  the  wall 
where  the  steepest  gradients  exist.  This  result  has  led  to  the  conclusion  that  the 
local  gradients  across  the  sensor  arrays  were  a  significant  cause  of  the  discrepancies 
between  the  x-type  probe  mean  data  and  the  "true"  profiles. 

Distribtuions  of  the  various  measurements  of  the  streamwise  and  vertical  velocity 
fluctuations  and  the  Reynolds  stress  distributions  are  given  in  Figs.  A-3  and  A-4. 

No  "true"  or  best  distributions  of  these  quantities  are  known  for  these  profiles  and 
so  the  data  at  a  given  location  in  a  profile  were  compared  with  the  average  of  all  the 
like  data  taken  at  that  location.  The  plots  of  Figs.  A-3  and  A-4,  then,  show  distri¬ 
butions  of  the  inconsistencies  between  the  various  measurement  techniques  for  the 
various  profiles.  The  agreement  between  the  separate  measurements  for  u'/Ue  and  v'/Ue 
were  generally  very  good  (±  *  1/2  percent  for  u'/Ue  and  ±  =  1/5  percent  for  v'/Ue). 

If  the  average  levels  of  u'/Ue  and  v'/Ue  are  approximated  as  0.06  and  0.04  respectively 
the  above  inconsistencies  are  equivalent  to  ±  8  percent  and  ±  5  percent  uncertainties 
in  the  fluctuating  velocities  (u1  and  v')  themselves.  Note  that  for  the  u'/lL  plot 
of  Fig.  A-3  data  were  included  for  a  "single  horizontal  wire".  These  particular 
data  were  obtained  with  an  analog  data  system  consisting  of  a  polynomial  linearizer 
and  a  true  RMS  voltmeter.  All  the  other  data  of  Fig.  A-3  were  obtained  with  the  pre¬ 
viously  described  analog-digital  system.  The  consistency  between  the  measurements 
for  these  two  very  different  techniques  is  excellent.  The  largest  inconsistencies  were 
observed  for  the  Reynolds  stress  distributions  of  Fig.  A-4.  This  is  not  at  all  sur¬ 
prising  as  it  is  far  more  difficult  to  measure  correlated  than  single  fluctuating 
quantities.  For  these  Reynolds  stress  data  a  ±  15  percent  inconsistency  band  encom¬ 
passes  nearly  all  the  measurements. 

Plots  of  the  average  "bias"  and  "standard  error"  for  the  fluctuating  quantities 
of  Fig.  I 1-3  and  I 1-4  are  presented  in  Fig.  A-5.  The  overall  profile  errors  for  the 
stress  measurements  were  clearly  much  larger  than  the  errors  for  the  individual  velo¬ 
city  fluctuations.  Unlike  the  mean  profile  results  of  Fig.  A-2  the  uncertainties 
of  the  fluctuating  quantities  for  the  various  boundary  layer  thicknesses  were  nearly 
constant. 
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Figure  A-3.  Discrepancies  Between  Various  Local  Measurements 
of  Fluctuating  Velocity  Components 
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Figure  A-4.  Discrepancies  Between  Various  Local  Turbulent 
Shear  Stress  Measurements 
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Figure  A-5.  Overall  Discrepancies  Between  Different  Fluctuating 
Measurement  Techniques  for  Various  Profiles 
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APPENDIX  B  -  EXPERIMENTAL  PROFILE  DATA 


All  data  for  each  profile  are  grouped  together  as  follows: 

a)  mean  profile  data 

b)  mean  profile  data  tabulation 

c)  boundary  layer  property  tabulation 

d)  fluctuating  profile  plots  (A-E) 

e)  fluctuating  profile  data  tabulation  (A-B). 

The  profile  data  are  presented  in  the  following  order: 
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.867 

.ei9 

-3  .266 

21.376 

16.282 

522. 5b7 

Hi 

.2853 

.972 

87.00 

78.16 

.679 

.615 

-2.971 

21.672 

16.158 

555.852 

92 

.3039 

.502 

88.15 

78.10 

.691 

.818 

-2  ,b88 

21.960 

18.275 

551.15? 

93 

.3163 

.526 

89.00 

73.57 

.900 

.651 

-2.873 

22.171 

19.011 

620. 168 

99 

.3332 

.551 

89.90 

78.06 

.909 

.821 

-2  .287 

22.397 

18.383 

685.213 

95 

.396  1 

.575 

90.60 

73.36 

.916 

.  868 

-2.028 

22.615 

19.297 

678.282 

96 

.3631 

.bTO 

91.59 

73.85 

•  9  2o 

.656 

-1.828 

22.61b 

19.170 

7C7.8fct 

97 

.3766 

•  626 

92.  38 

73.03 

.938 

.888 

-1.629 

23.018 

19.786 

738.28? 

96 

.3932 

.650 

93.29 

72,99 

.983 

.887 

-1  .815 

23.226 

19.600 

766.10* 

99 

.9065 

.675 

93.96 

73.00 

•  9  5C 

.886 

-1  .231 

23.812 

15.757 

755.816 

50 

.9233 

.700 

99.66 

72.62 

.957 

.909 

-1  .062 

23.562 

20.312 

828.750 

51 

.9369 

.725 

95.16 

72.36 

.962 

.926 

-.938 

23.706 

20.671 

859.168 

52 

.9537 

.750 

95.62 

72.11 

.969 

.981 

-.772 

23.672 

21.012 

6  6  3  •  c  7  6 

53 

.9669 

.779 

96.27 

72.28 

.973 

.933 

-  .660 

23.583 

20.638 

912.616 

59 

.9833 

.799 

96.77 

72.  36 

.978 

.925 

-.536 

28.108 

20.666 

981 . 688 

55 

.9966 

.829 

97.32 

72.00 

.988 

.988 

-.396 

28  ,28b 

21.163 

971.853 

5b 

.5133 

.898 

97.61 

71.83 

.967 

.956 

-.326 

28.317 

21.399 

18§8:K? 

57 

.5287 

.879 

97.97 

71.57 

.990 

.978 

-.237 

28.807 

21.752 

58 

.5933 

.898 

98.22 

71.78 

.993 

.  968 

-.176 

28.867 

21.527 

1056.539 

59 

:*»H 

.5886 

.923 

98.96 

71.67 

.995 

.968 

-.118 

28.530 

21.624 

1087.763 

60 

98.67 

71.55 

.997 

.975 

-.062 

28.581 

21.786 

1116.566 

61 

96.77 

71.88 

.996 

.980 

-.038 

28.606 

21.882 

1146.798 

62 

.6239 

.997 

98,77 

71.05 

.996 

1.006 

-.038 

28.605 

22.872 

1 1 75.628 

63 

.6183 

1.022 

98.86 

71.81 

1 .000 

.968 

-.01.0 

28.638 

21.976 

1208.657 

69 

.6339 

1.C97 

96.92 

71.87 

l.OCG 

.980 

-.002 

28.682 

21.693 

1238.076 

65 

.6968 

1.072 

98.97 

71.98 

1.000 

.960 

.011 

.027 

28.655 

21.686 

1268.075 

66 

.6639 

1.097 

99.03 

71.58 

1.0C1 

.973 

28.670 

21.735 

1253.857 

67 

.6789 

1.121 

99.07 

71.25 

1.0C1 

.998 

.036 

28.680 

22.197 

22.264 

1321.787 

68 

.6933 

1.196 

99.02 

71.20 

1.001 

.997 

.028 

28.666 

1350. 776 

69 

.7083 

1.171 

98.96 

71.50 

l.CCC 

1.009 

.006 

28.652 

22.580 

1379.959 

70 

.7231 

1.196 

99. ce 

71.33 

1.002 

.989 

.038 

28.662 

22.090 

141C.  157 

71 

.7386 

1.221 

99.03 

70.91 

1.001 

1.015 

.026 

28.670 

22.668 

1839.821 

72 

.75  39 

1.295 

96.99 

70.88 

1.001 

1.017 

.998 

.016 

28.659 

22.707 

1867.665 

73 

.7687 

1.271 

99.02 

71.18 

1.001 

.028 

28.666 

22.266 

1897.673 

79 

.7839 

1.295 

96.99 

71.21 

l.COO 

.996 

.005 

26.685 

22.252 

1526.312 

75 

.7989 

1.320 

99.05 

71.19 

71.52 

1.001 

.997 

.031 

26.675 

22.273 

1555.536 

76 

.8135 

1.395 

99.11 

1.CC2 

.977 

.088 

26.651 

21.628 

1588.958 

77 

.1265 

1.369 

98.99 

70.86 

1.001 

1.018 

:8  hi 

26.661 

22.727 

1618.176 

78 

.89  33 

1.399 

98.95 

71.17 

1.000 

.998 

26.650 

22.299 

1683.012 

79 

.8567 

1.919 

96.96 

71.23 

l.COC 

.995 

.009 

26.652 

26.657 

firm 

1673.015 

1C 

.1733 

1.993 

98.96 

70.89 

1.001 

1.016 

.018 

1701.855 

81 

.8887 

1.969 

99.09 

71.52 

l.gci 

.977 

.031 

24.678 

21.630 

1731.862 

6§ 

i-Mll 

1.999 

98.99 

71.96 

1.0C1 

.981 

:Sos 

28.661 

21.914 

1760.886 

2712.738 

i:i\l 

IV.lt 

71.29 

1.000 

.992 

24.652 

22.168 

8  8 

1.9533 

71.21 

1.000 

.997 

.009 

28.653 

22.256 

3605. 562 

85 

2.9763 

9. 096 

96.91 

71.51 

l.OCQ 

.578 

-.003 

26.681 

21.636 

4826.350 

86 

3. 0033 

9.969 

99.01 

71.69 

1.001 

.970 

.021 

26.665 

21.657 

5851.216 

Table  B1 


RUN  NO.  3.  POINT  3. 


BOLNDARY  LAYER  PROPERTIES 

standard 

LINEAR 

SUBLAYER 

interpolation 

FUNCTION  F 

TO  WALL 

WALL  TO  Y ♦ 

FREE  STREAM  VELOCITY 

= 

98.922 

98.922 

FREE  STREAM  TEMPERATURE 

z 

71 . 1 5C 

WALL  TEMPERATURE 

z 

87.410 

WALL  HEAT  FLUX 

• 

.05120 

FREE  STREAM  OENSITY 

z 

.07523 

FREE  STREAM  KINEMATIC  VISCOSITY 

z 

.C00162e 

DENSITY  OF  FLUID  AT  WALL 

z 

.07300 

KINEMATIC  VISCOSITY  OF  FLUID  AT  WALL 

z 

•  COO  1 7  1 7 

WALL/FREE  STREAM  DENSITY  RATIO 

z 

.97028 

LOCATION  REYNOLDS  NUMBER  «  R  E  X  » 

z 

2633499.75 

INPUT  VALUE  OF  VELOCITY  DELTA 

z 

.61000 

INPLT  VALUE  OF  TEMPERATURE  DELTA 

z 

.67000 

CALCULATED  delta 

z 

.60504 

DELTA  99. SX  INPUT 

z 

.00000 

DISPLACEMENT  THICKNESS  IDELSTAR) 

z 

.09026 

.09063 

MOMENTUM  THICKNESS  (THETA) 

z 

.06280 

.06260 

ENERGY-DISSIPATION  THICKNESS 

z 

.11090 

.11076 

enthalpy  thickness 

z 

. D02  74 

.00273 

SHAPE  FACTOR  12  (DELSTAR/THETA) 

z 

1.43730 

1.44315 

SPAPL  FACTOP  32  (ENERGY/THETAI 

z 

1.76594 

1.7641G 

MOMENTUM  THICKNESS  REYNOLDS  NUMBE  P 

z 

3180.47 

3160.34 

DISPLACEMENT  THICKNESS  PEYNOLOS  NUMBER 

z 

4571.29 

4589.72 

SKIN  FRICTION  COEFFICIENT 

z 

.003195 

FRICTION  VELOCITY 

z 

4.01412 

LAW  OF  THE  WALL  CONSTANT  (K) 

z 

. 4 10DC 

LA*  OF  THE  WALL  CONSTANT  (CJ 

z 

5.0.00  00 

WAKE  STRENGTH 

1 

.49062 

CLAUSEPS  ’DELTA*  INTEGRAL 

Z 

-2.07555 

-2. 16672 

CLAUSERS  *  G ’  INTEGRAL 

z 

14.6892 2 

14.91446 

DISPLACEMENT  THICK Nt.SS  -  CONSTANT  DENSITY 

z 

.0e589 

.08792 

MOMENTUM  THICKNESS  -  CONSTANT  DENSITY 

z 

.06337 

.06336 

SHAPE  FACTOR  12  -  CONSTANT  DENSITY 

z 

1.35544 

1 .3675e 

LOCATION  -X' 

m 

52.00000 

Te  =  0.2% 
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fTt  ~*<r- 


IX  II 


-  uV/uV' 


Figure  B2A.  Boundary  Layer  Turbulence  Quantities 


UV  TRIPLE  PRODUCTS 


Figure  B2C.  Boundary  Layer  Triple  Product  Distributions 


B2D.  Boundary  Layer  Skewness  and  Triple  Product  Correlation 
Coefficient  Distributions  x  =  52  in,  Te  =0.2% 


Fluctuating  Profile  Data 


x  =  52  in,  Te  =  0.2% 


r  : 

7/ 
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N 
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3 
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it 
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.27  6 

.0567 

.0362 

.0923 
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.060*. 

c 

•  2CC 

.370 

.0527 

.0357 

.09)9 

.0351 

.  79  79 

.5592 

.  C  5  9 

6 

.250 

.96  3 

.0976 

.0327 

.0377 

.0312 

.SSol 

.5996 

.ct  e  = 

7 

.3  CO 

.556 

.0909 

.0269 

.0329 

.0272 

.9965 

.  5951 

•  etc: 

6 

.350 

.69  E 

.0360 

.0251 

.029$ 

.0225 

.  30b8 

.5171 

.  05  2  9 

9 

.9  CO 

.791 

.0291 

.0210 

.0239 

.0179 

.1837 

.9712 

.09  31 

1C 

.950 

.  e  3  3 

a  C2  1  9 

.0165 

.0189 

.0129 

.  0931 

.9525 

.03  2' 

11 

.SCO 

.926 

.0191 

.0126 

.0123 

.0060 

.0391 

.  3967 

•  022‘ 

12 

.5  50 

1.015 

.0061 

.0090 

.0079 

.0095 

.0123 

.2911 

.01  “9 

13 

.6  CO 

1.111 

.0055 

.0069 

.  0098 

.0026 

.0090 

.1959 

.009’ 

1<* 

.650 

1.209 

.0092 

.0050 

.0035 

.0016 

.0015 

.1296 

.  00  7- 

IS 

•  7  CO 

1.296 

.0037 

.0038 

.002$ 

.0011 

.0007 

.0627 

.0059 

16 

.750 

1.389 

.0C3! 

.0032 

.0021 

.0009 

.0005 

.0733 

.00  50 

17 

l.OCO 

1.852 

.0028 

.0021 

.0019 

.0009 

*.0001 

-.0216 

.00  36 

is 

3.CC0 

5.556 

.0022 

.0016 

.0011 

.0005 

-.0001 

-.0996 

.  00  3: 

y  : 

Y/ 

Y  T  A  U  / 

VC/ 

V . T ■ PNC* 

K 

INCHES 

OELTA 
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6’2/C 2 

UTAO 

UT  Ac. 

C  .  6  w  L 

3.CC0 

.050 
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.2195 

.3261 
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•AUi 

.0019 

3 
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.2168 
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9 
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.9971 
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.  3527 
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5 

.  2  CO 

.  3  7  C 

.2192 
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.2215 
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.2376 

6 

.2  50 

.9b3 

.20s3 
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.2253 

.2955 

.  5053 

.2895 

.6056 

7 
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.556 

.2058 

.9662 

.  2328 

.3009 

.6197 

,39eG 

.  756  = 

8 

.350 

.69  8 

.1808 

.9  637 

.2250 

.3113 

.  7033 

.3589 

•  56  So 

9 

.900 

.79  1 

.1629 

.9558 

•.2369 

.3079 

.8*36 

.9027 

.961- 

1C 

.950 

.833 

.1935 

•  929e 

.2539 

.3163 

.9265 

.9266 

.  36  6  = 

11 

.500 

.926 

.1266 

.3919 

.3123 

.2962 

.8932 

.3730 

•  2328 

12 

.550 

1.019 

.0971 

.3195 

.3es5 

.2960 

.6531 

.2312 

.192- 

13 

•  6CC 

1.111 

.0693 

.3225 

.9307 

.2968 

.  3620 

.1171 

.  0756 

19 

.650 

1.209 

.0966 

.329? 

.  99  77 

.2280 

.1237 

.0915 

.00  26 

IS 

.7  CO 

1.296 

.0355 

.3902 

.9273 

.1625 

-.0559 

,0229 

.007  = 

16 

.750 

1.389 

.0331 

.9299 

.9095 

.1706 

-.0993 

.0115 

.0021 

17 

l.OCO 

1.852 

-.0109 

.5391 

.3193 

.1916 

-.0131 

.0007 

-.000- 

18 

3.0C0 

5.556 

-.0255 

.5770 

.2831 

.1399 

.0590 

.0006 

-.cco? 

N 

T  S 

INCHES 

Y/ 

DELTA 

V  ’T •/ V ’T  • 

T • /TTAU 

T’/ITk-TEI 

*16 

G 1 6 

PR  T 

\ 

3:8S8 

5.556 
.09  3 

•:H!3 

2:m 

:8?li 

-.8704 

.8526 

Ml 

-1.5202 

3 

•  ICO 

.165 

.5389 

1  .6112 

.0714 

.7351 

.0312 

6 

.150 

.27  6 

.5909 

1  .9479 

.0641 

.7457 

.0327 

5 

.200 

»37C 

.5535 

1.3417 

.0594 

.6058 

.  0426 

6 

.963 

.5609 

1  .2262 

.0543 

.6372 

.0599 

.6707 

7 

.556 

.5726 

1.1399 

.0505 

.8978 

.0654 

.6062 

8 

.350 

.69  e 

.5697 

1  .02C9 

.  0452 

.9203 

.0600 

.6041 

9 

.900 

.79  1 

.5780 

.9421 

.0417 

1.0504 

.  1096 

.6456 

10 

Mi 

.5789 

.6292 

.0367 

1.2202 

1.4047 

.1562 

.5457 

11 

.5389 

.6800 

.0301 

.1718 

.4556 

12 

.550 

1.015 

.9797 

.5065 

.0224 

1  .6409 

.1820 

.3898 

13 

.600 

1.111 

.9029 

.3464 

.0153 

2.2637 

.1622 

.2666 

is 

.650 

.700 

1.209 

.2980 

.2164 

.0096 

1 .7649 

.1272 

.2274 

1.2ft 

.1955 

.1539 

.0068 

1.2742 

.0716 

.2876 

16 

.750 

1.366 

.0658 

.1294 

.0057 

.5113 

.0268 

.4156 

17 

l.OCO 

1.852 

-.0108 

.1200 

.1123 

.0053 

-.2405 

-.0017 

-1.6186 

16 

3.0C0 

5.55  6 

-.0075 

.0050 

-.0799 

.0000 

-3.0404 

Table  B3A 


Fluctuating  Profile  Data 


x  =  52  in.,  Te  =  0.2‘: 


Y  : 

Y  / 

C  *  24*/ 

U,6,2/ 

U  •  2  v  *  / 

V'V’2/ 

K 

iNCnLi 

DEL  T  A 

0  '.J/UC  3 

UE  3 

* ' 3/CC 3 

UE  3 

yC  3 

UE  3 

V’  3/1  E  : 

X  1CCCC 

X  10000 

x  10PD0 

X  10000 

X  lOOCO 

X  100CC 

X  1CCC' 

* 

l:m 

-:8£c2 

:?£££ 

-•m 

:8£B£ 

-:sm 

.  CC  00 
.  15  35 

3 

.  1  CO 

.16  5 

-.4256 

-.0768 

.0236 

.0188 

.2068 

-.1456 

.05  35 

4 

.  1  so 

.27  6 

-.4271 

-.0677 

.0730 

-.1162 

.2195 

-.1677 

.  12t5 

5 

•  2  CO 

.  3  7  C 

-.5037 

-  .0664 

.0291 

-.1354 

.2366 

-.1623 

.116  5 

6 

.2  50 

•  4  b  7 

-.4229 

-.0693 

-.0218 

- .1587 

.1990 

-.1604 

.14  33 

7 

.  3CC 

.556 

-  .4237 

.0058 

.0337 

-.1068 

.1947 

-.1525 

.176- 

8 

.  350 

.64  £ 

-.3132 

-.0047 

•  004  n 

-.1094 

.1453 

-.1200 

.05  5? 

9 

.4  00 

.74  1 

-  .2465 

.0198 

.0222 

-.1  IOC 

.1064 

-.0873 

.076'’ 

10 

•  4  SC 

.83  3 

-  .  1 2C 4 

-.0067 

-.0138 

-.0653 

.0566 

- .  C  5  C  7 

.  C  4  7  ** 

11 

.SCO 

.926 

-.0407 

-  .0060 

-.0038 

-.0185 

.0211 

- . C  2  06 

.0225 

12 

•  s  so 

1 .01  9 

-  .0077 

.0022 

.0014 

-.0057 

.0039 

-.0049 

.  CC6- 

13 

.6  00 

1.111 

-.0013 

-.0000 

.0001 

- .0006 

.0009 

-.0009 

.00  21 

14 

.  .650 

1.204 

.0001 

-.0001 

-.0002 

.0001 

.0002 

-.0001 

.  OC  1  7 

IS 

.  7  CO 

1.296 

.0001 

- .0000 

-.0000 

.0001 

.  oocn 

.0000 

.  00  2' 

16 

.  7  SC 

1.385 

•  COoC 

-  .0000 

.OCCO 

.COCO 

•  u  0  Ci  C 

•  COOC 

.  0  C  C  1 

17 

1.000 

1.652 

-.oc5o 

.crcc 

.0000 

-.0000 

.0000 

.COCO 

.ccc; 

18 

3.000 

S.SS6 

-  .COCO 

- .0000 

.COCO 

-.0000 

.  OOuC 

.coco 

.  0  C  0  0 

N 

INCHES 

oel{* 

SU  V 

SU  H 

SV 

Sw 

ST 

PSIEV 

PSIE. 

1 

3.  0  CO 

5.556 

- .0085 

.0000 

.1593 

.0000 

-.3002 

-.0024 

.  C  C  u  ' 

2 

.0  50 

.09  3 

-.1185 

.0052 

.2349 

.1338 

.0078 

.1497 

-.015' 

3 

.ICC 

.165 

-.2137 

.0304 

.1462 

.0263 

-.0051 

.1506 

-.  0  3  8  4 

4 

.150 

.278 

-.2698 

-.2203 

.2255 

.0999 

.052e 

.1961 

-.0444 

.  2  CC 

.  3  7  C 

-.4216 

-.2672 

.2600 

.0429 

.0628 

.2753 

-.05  le 

6 

.250 

.463 

-.5135 

-.5573 

.4093 

-.0416 

.2129 

.3215 

-.0661 

7 

.300 

.556 

-.7500 

-.4706 

.5202 

.0915 

.2332 

.4568 

.0103 

8 

.3*o 

•64  e 

-.9242 

-.7918 

.6275 

.0063 

.3772 

.5525 

-.020' 

9 

.  4  CC 

.74  1 

-1  .392? 

-1.0619 

.6290 

.1506 

.  5875 

.7468 

.0767 

10 

.450 

.833 

-2.0364 

-1  .6232 

1.0497 

-.2327 

.9722 

1.0512 

-.0E3’’ 

11 

.SCO 

.92  6 

-2.6934 

-2.5342 

1.139P 

-.2045 

1.3906 

1.2750 

-.169? 

12 
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VELOCITY  AND  TEMPERATURE  RATIOS 


VELOCITY  AND  TEMPERATURE  DISTRIBUTIONS  IN  UNIVERSAL  COORDINATES 


7S— 12— 100- 1 


Mean  Profile  Data 
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RUN  NO.  3.  POINT  2. 

BOUNDARY  LAYER  PROPERTIES  STANDARD 

LINEAR  SUBLAYER 

INTERPOLATION  FUNCTION  FRO* 

TO  'WALL  •> ALL  TO  Y*  =  25 

FREE  STREAM  VELOCITY  r  99.17i  99.171 

FREE  STREAM  TEMPERATURE  =  73.270 

WALL  TEMPERATURE  =  90.27C 

WALL  HEAT  FLUX  :  .05020 

FREE  STREAM  OENSITY  =  .07493 

FREE  STREAM  KINEMATIC  VISCOSITY  =  .0001629 

OENSITY  OF  FLUID  AT  »ALL  =  .07262 

KINEMATIC  VISCOSITY  OF  FLUID  AT  WALL  =  .C001733 

WALL/FREE  STREAM  OENSITY  RATIO  =  .96909 

LOCATION  REYNOLDS  NUMBER  (REX)  =  3426173.22 

INPUT  VALUE  OF  VELOCITY  DELTA  =  .90000 

INPUT  VALUE  OF  TEMPERATURE  DELTA  =  .97000 

CALCULATED  DELTA  =  .61093 

DELTA  99.5*  INPUT  =  .OOOCP 

displacement  thickness  (delstar)  =  .12319  .123C6 

momentum  thickness  (theta)  r  .08551  .08564 

ENERGY-DISSIPATION  THICKNESS  =  .15118  .15l4e 

ENTHALPY  THICKNESS  =  .00370  .00371 

SHAPE  FACTOR  1?  ( 0 EL S T AR / T H E T A  )  :  1.44071  1.43351 

SPAPE  FACTOP  32  (ENERGY/THETA  I  =  1.76804  1.76469 

MOMENTUM  THICKNESS  REYNOLDS  NUMBEP  =  4310.75  4327.66 

DISPLACEMENT  THICKNESS  REYNOLDS  NUMBER  =  6210.57  62D3.78 

SKIN  FRICTION  COEFFICIENT  =  .0D2940 

FRICTION  VELOCITY  =  3.66252 

LA-  OF  THE  WALL  CONSTANT  (K>  =  .41000 

LA*  OF  The  WALL  CONSTANT  (C)  =  5.0H000 

WAKE  STRENGTH  =  .57972 

CLAUSEPS  ’DELTA •  INTEGRAL  =  -2.88141  -3.06592 

CLAUSERS  ’G’  INTEGRAL  =  21.95430  21.63072 

DISPLACEMENT  THICKNESS  -  CONSTANT  DENSITY  =  .11569  .11941 

MOMENTUM  THICKNESS  -  CONSTANT  DENSITY  =  .08625  .08660 

SHAPE  FACTOR  12  -  CONSTANT  DENSITY  =  1.34363  1.37890 


LOCATION  -X-  6e. 00000 

Te  =  0.2% 
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An/An-O]  \  V  L  en/,A  □ 


UV  TRIPLE  PRODUCTS 


Boundary  Layer  Triple  Product 


Boundary  Layer  Skewness  and  Triple  Product 
Coefficient  Distributions  x  =  68  in  T  = 


Figure  B4E.  Boundary  Layer  Flatness  Distributions 


Fluctuating  Profile  Data 
x  *  68  in.,  Te  =  0.2% 
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Fluctuating  Profile  Data 
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FREE  STREAM  TEMPERATURE 
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Table  B8 


q/q^  MEAN  PROFILE 


Figure  B6B.  Boundary  Layer  Turbulence  Quantities, 
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-.9180 

-1.1102 

-.8277 

-.3179 

.0992 

.0656 

-.0162 
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*:8I5S 

.0653 
.0373 
.0955 
.1603 
.1  397 
.3637 
.3005 
.6737 
.7085 
.8167 
.8767 
.3672 
.0386 
.0599 
.0339 
.0000 


_.»«* 
2.7397 
2.8000 
2.8636 
3. 07b6 
3.1711 
3.2803 
3.5625 
5.0039 
8.5«60 

H-.S1S1  1) 

12.6789 
3.7858 
3.3259 
2.9652 

•  ooco 


3.2235 
2.6698 
2.6266 
2.8050 
2.7106 
2.8666 
2.8528 
3.1062 
3.3716 
6.2738 
6.6175 
.6509 
.2666 
11.5225 
16.3031 
3.0096 
2.8138 
.0000 


I:  ill!  i:!88S 


3.1962 
3.0816 
3.1670 
3.2216 
3.2629 
3.6158 
3.9022 
6.8173 
5.8260 
9.5023 
12.2509 
10.9790 
6.2007 
3.1608 
3. 0069 
.0000 


_ 09b 

3.0630 

3.9321 

3.1262 

3.1896 

5.1661 

3.5636 
3. 5012 
6.1965 
5.6*09 
10.0596 
16.5*81 
15.6552 
6.6586 
2.9333 
2.6356 
.CCOC 


2.675' 

9 ,9t  ;~ 

3. 702c 
j.67c: 
3.952’ 
3.511' 
3.9390 

3.3’:^ 
3.995* 
3.739* 
9.702? 
8.999? 
19.559* 
3  2 . 5  *  1 5 
92.5511 
2.892* 
2 . 6 2  06 
.coo: 


Table  B9B 


,<*  v  — 


Mean  Profile  Data 


x  =  52  in.,  Te  =  1 . 82 


Y 

N  X KCHc  5 

1  .0063 

2  .0095 

3  .010“ 

m  .one 

5  .0125 

6  .0142 

7  .0157 

8  . 01 b  7 

9  .0167 

10  .0207 

11  .0224 

12  .0241 

13  .025b 

14  ,0274 

15  .0297 

16  .0315 

17  .0332 

18  .0*96 

19  .04  o  3 

20  .0537 

21  .0596 

22  .0666 

23  .0737 

24  .0795 

25  .0866 

26  .0936 

27  .0996 

28  .1065 

29  .1135 

30  .1193 

n  :im 

S3  .15^6 

34  .1662 

35  .1855 

36  .2036 

37  .2203 

38  .2386 

39  .2554 

40  .2733 

41  .2904 

42  .3085 

43  .3535 

44  .3984 

45  .4435 

46  .4884 

47  .5333 

48  .5784 

49  .6234 

50  .6684 

51  .7133 

52  .7583 

53  .8034 

54  .8484 

55  .*’35 

56  .93o4 

57  .9834 

58  1.0283 

59  1.0734 

60  1.1184 

61  1.1635 

62  1 .20o4 

63  1.9362 

64  2.6084 

65  3.3087 

66  4.0084 


Y/ 

OElT  4 
.010 
.011 
.512 
.314 
.014 
.016 
.016 
.019 
.0  22 
.024 
«S2o 
.026 
.030 
.0*2 
•  03  4 
.036 
.036 
.046 
.053 
.062 
.069 
.0*7 

.ces 

.092 

.ICC 

.108 

.116 

.123 

.131 

.138 

.146 

.154 

.174 

.194 

.214 

.235 

.254 

.275 

.255 

.315 

.335 

.356 

.406 

.459 

.511 

.563 

.615 

.667 

.719 

.771 

.823 

.874 

.926 
.978 
1.030 
1.562 
1.134 
i  .let 
1.238 
1.290 

1.342 

1.393 

2.200 

3.008 

3.815 

4.622 


u 

FT/S;.C 
45.17 
4  7. 06 
4  8.69 
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53.  C7 
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55. 32 
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57.  Oo 
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68.22 
68.72 

69. 44 

70. 20 
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71.44 
71.67 

73.49 
74. 34 

75.82 

76.75 
77. 6b 

79.  C3 

80.  CO 

61.09 

81.56 

62. 82 
85.04 

87.00 

88.78 

90.26 
91.61 
93. C3 

94. 22 

95.27 
95.90 
96.55 

97.19 
97. 39 
97.60 
97. 89 
98.01 
97.97 
98.  15 

98.19 

98.22 
98.19 
98.25 
98.12 
9  8  .  1 0 
98.01 


T 

occ.r 

8C.92 

34.89 
84.47 
84. 04 

87.78 
S3.  12 
82.66 
8?. 58 
82.67 
62.23 
61.45 
81.62 

81.40 
81.09 

80.69 
8  0.82 
sc. eo 

80.89 
80.51 

79.40 
79.54 
78.76 
79. C3 
79.04 

78.10 
78.72 

76.42 
76.31 
77.80 
77.66 
7«.0t 
77.87 
77.05 

76.8  4 

76.70 

76.10 

75.62 

75.8  0 
76.10 
75.64 

75.42 

75.16 

74.74 

73.75 
77.57 

72.71 

72.69 
72.24 

72.42 

72.35 

71.69 
71.28 

71.63 

71.36 

71.37 
71. C2 
71.14 
70.83 
71.19 
71.61 
71.13 
70.68 
71.41 

71.16 
71.33 
71.22 


U/UC 
.460 
.480 
.4  96 
.507 
.515 
.523 
.541 
.542 
.556 

.564 

.570 
.577 
.582 
.569 
.592 
•  6  0  0 
.602 
.615 
.629 
.640 
.653 
.661 
.672 
.679 
.668 
.695 
.  7  Cl 
.708 
.716 
.721 
.728 
.731 
.749 
.756 
.773 
.762 
.794 
.806 
.816 
.627 
.831 

.844 

.867 
.687 
.905 
.920 
.934 
.946 
.960 
.971 
.978 
.984 
.991 
.993 
.995 
.996 
.999 
.999 
1.000 
1.QC1 
1.001 
1.001 
1.002 
1.000 
1  .ooc 
.999 


ThET4 

.377 

.420 

.437 

.455 

.466 

.494 

.513 

.516 

.513 

.531 

.563 

.556 

.565 

.578 

.587 

.590 

.590 

.586 

.602 

.  64  B 

.643 

.675 

.664 

•  664 

.703 

.677 

.689 

.  b94 

.715 

.712 

.704 

.712 

.747 

.755 

.761 

.786 

.606 

.796 

.786 

.805 

.614 

.  824 

.843 

.864 

.891 

.927 

.920 

.947 

.939 

.942 

.970 

.987 

.972 

.983 

.983 

.997 

.993 

1.006 

.990 

.973 

.993 

1.012 

.961 

.992 

.915 

.989 


U-UE 

-13.631 
-13.193 
-12.726 
-12.449 
-12.263 
-11.892 
-11.598 
-11.564 
-11.225 
-11.016 
-10.858 
-10.695 
-10.563 
-10.384 
-10.306 
-10.115 
-10.058 
-9.731 
-9.364 
-9  .083 
-8.7  60 
-6 .566 
-6  .261 
-8  .1 16 
-7 .884 
-7 .694 
-7 .566 
-7.362 
-7.186 
-7  »  D4  4 
-6 .865 
-6.607 
-b .337 
-6.119 
-5.739 
-S  .496 
-5.212 
-4  .911 
- 4 .661 
- 4 .361 
-4  .260 
-3.935 
-3.363 
-2.859 
-2 .401 
-2 .019 
-1.670 
-1.306 
-1 .ODD 
-.729 
-.566 
-.400 
-.235 
-.183 
-.128 
-.054 
-.024 
-.034 
.013 
.022 
.030 

.024 

.039 

.006 

.000 

-.024 


U  (  *  ) 
11.632 
12.120 
12.536 
12.615 
13.0C1 
13.371 
13.866 
13.700 
14.035 
14.247 
14.406 

14.569 
14.70C 
14.000 
14.956 
15.149 
15.206 
15.532 
15.899 
lb. 180 
16.504 
1 6 • o9  6 
lb. 96  3 
17.147 
17.375 

17.570 
17.656 

17.682 
18.076 
16.220 
16.395 
18.457 
18.927 
19.14s 
19.525 
19.765 
20.052 
20.353 
2  0  .  bO  3 

20.683 
21.004 
21.326 
21.901 
22.405 

22.863 
23.244 
23.593 
23.956 
24.764 
24.534 
24.698 

24.864 
25.025 
25.081 
25.136 
25.209 
25.240 
25.229 
25.276 
26.286 
25.294 
25.287 
25.302 
25.269 
25.264 
25.239 


T  (  ♦  ) 

7 .706 
6.576 
6.932 

9.297 
9.524 
1C. 087 
1C. 471 
10.542 
10.466 
10.635 
11.507 
11.363 
11.S4E 
1  1.606 
11.562 
12 .042 
12.054 
11.576 
12.30C 
13.2«2 
13.124 
13.793 
13.559 
13.555 
14.354 
13.622 
14.075 
14.170 
14.603 

14.536 
14 .3eu 
14.546 
15.246 
15.421 
16.542 
lb. 046 
lb. 457 
16.305 
16.056 
16.442 
16.627 
lb. 835 
17.2D7 
16.051 
18.205 

16.536 
18.7e2 

15.334 

19.183 

19.241 

19.799 

20.149 

19.653 

20.079 

20.073 

20.370 

20.269 

20.536 
20.222 
19.867 
20.275 
20.657 
20.037 
20.254 
20.106 
20.199 


Y  I  •  ) 
16.35! 
17.565 
15.224 
2 1  •  6  C  4 
2 3 • C4  c 
26.226 
26.55* 
30.636 
34.522 
36.7Ce 
41.342 
44.475 
47 . ?4f 

50.556 
54.757 
56.114 
61. ?4? 

73.544 

85.353 
95.032 
1  09 . 5  C  * 
122.605 
136. e55 
146.566 
16C.C41 
172.574 
183.633 
196.351 
209. 253 
215. 5-3 
233.356 
246.116 
077.63! 
310.073 
341.555 
375.320 
406.101 
4J5. P3C 
470.55 
503.765 
535.305 
5b6.ttfc 
651 .frCS 
734.365 
817.451 
9  CC • ?4  5 
983. CC6 
106b. 132 
1149. 074 
1232.015 
1314.773 
1597.714 
1480.640 

1563.762 
1646.906 
1729.665 
1812. 807 
1895. 364 
1978.490 
2061.43? 
2144.556 
2227. J 1 5 
3517.150 
48C7. 722 
6098.479 

738e. 13C 


Table  BIO 


{jpg**  -v 


RJN  NO •  3.  POINT  3. 

BOlNnARY  layer  properties  STANDARD 

linear  Sublayer 

INTERPOLATION  FUNCTION  PRC*' 

TO  VALL  -ALL  TO  Y*  =  I5 


FPEr  STREAM  VELOCITY 

z 

70.96& 

96.101 

FREE  STREAM  TEMPERATURE 

z 

WALL  TEMPERATURE 

z 

9*)  .980 

WALL  HEAT  FLUX 

z 

.07908 

FREE  STREAM  DENSITY 

z 

.07541 

FREE  STRFAM  KINEMATIC  VISCOSITY 

z 

. COO  162  3 

DENSITY  OR  FLUID  AT  WALL 

z 

.07214 

kinematic  viscosity  of  fluid  at  WALL 

z 

. C0C175E 

•ALL/FREE  STREAM  DENSITY  RATIO 

z 

. 95669 

LOCATION  REYNOLDS  NUMBER  (REX) 

z 

2618495.47 

INPUT  VALUE  OF  VELOCITY  DELTA 

z 

1.00000 

INPUT  VALUE  OF  TEMPERATURE  DELTA 

z 

1 .00000 

CALCULATED  DELTA 

z 

.86723 

DELTA  95.5*  INPUT 

z 

.00300 

DISPLACEMENT  THICKNESS  (UELSTAR) 

z 

. 124  06 

.12411 

MOMENTUM  THICKNESS  (THETA) 

z 

.08716 

.08756 

ENERS Y-DISSIP ATION  THICKNESS 

z 

.  15548 

.15561 

ENTHALPY  THICKNESS 

z 

.00436 

.00437 

SHAPE  FACTOR  1?  (OELSTAR/THETA  ) 

z 

1.42349 

1.41709 

SHAPE  FACTOR  3 2  < ENERGY /7H£T A > 

z 

1.78378 

1  .  7  7  9  C  7 

momentum  thickness  Reynolds  number 

z 

4389,16 

4410.12 

displacement  thickness  Reynolds  number 

z 

6247.94 

6249.55 

SKIN  FRICTION  COEFFICIENT 

z 

.002998 

FRICTION  VELOCITY 

z 

3.86308 

L  A ■  OF  THE  WALL  CONSTANT  (K) 

z 

.41000 

la.  OF  The  WALL  CONSTANT  <C> 

z 

E.OODCD 

wAkE  strength 

- 

.46566 

CLAUSERS  ’DELTA*  INTEGRAL 

z 

-2.81396 

-3.D25G1 

CLAUSERS  INTEGRAL 

z 

2  2. 200  2  C 

19.93622 

DISPLACEMENT  THICKNESS  -  CONSTANT  DENSITY 

z 

.11555 

.11974 

MOMENTUM  THICKNESS  -  CONSTANT  DENSITY 

z 

.06806 

•  OS  8  S  D 

SHAPE  FACTOP  12  -  CONSTANT  DENSITY 

• 

1.31209 

1 .35299 

LOCATION  -X 

_ 

52.00000 

Te  =  1.8% 


Table  Bll 


-  uVAi'v' 


UV  TRIPLE  PRODUCTS 


Figure  B8C.  Boundary  Layer  Triple  Product  Distributions 


Figure  B8E.  Boundary  Layer  Flatness  Distributions  x  =52  in,  Te  =  1.8% 
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x  =  52  in. ,  Te  =  1 -8* 
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Fluctuating  Profile  Data 
x  =  52  in. ,  Te  =  1 .8" 
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■1 .Ox  ;x 

-.66x9 
-.  xS  7C 

-  .  1  7  ;9 
- .  Cxx  1 
-.0156 
-.ul7x 


■•0CX6 

■•otss 

:rifl 

.1124 

..ex? 
• .  06  9  S 
“  ,  OCEC 
.:  j'2 

..it? 
-.0100 
.  cr  3 1 


-1.2939 

•l.lStS 

-.7209 

-.3s  ;c 

::SSI? 

*.WwO 

-. is»: 
.ococ 


-.7771 

-  .  £9  r  P 

•i.::i: 

•l.x’15 

■1.7X21 

•  1 .OS  3  7 

.cccc 

.  c  c 

-  .  16  S  9 

.*occ 


x  .2  7  s 
.0=7 
.12$ 

.:c: 

.275 
.35  3 
.  x  *  6 

•  S  3S 

•  6x2 

.7u« 
.656 
.967 
i.:7; 
1.176 
1.263 
l  .51; 
x  .27* 

.:oc 


.1667 
.  3S  7* 
•  3S9E 
.507. 
«  7X  xl 

:«!« 
.  3771 
.1261 
.25x2 
.0X2! 
.Cxi* 
.0691 
.0000 


•.IS  1C 
•  07xS 
.0*02 
.  0021 

:$«* 
.  12  5X 
.3"39 
.  xC  10 
.3562 
.  3?6  0 
.2-9* 
.2736 

-:ixcit 

.ore: 

:&6cc 


.•  3/C-27 

loco: 


•  «J  U  i  X 

.  0  73  9 
-.2*76 

•  C 63  c 

.  .-6i, 

I  HCfe  r 
-.07X9 
.0*9  = 
.000* 


y  / 

Oi.  u 

Sc  V 

SC  H 

3  V 

X  .27  * 

.0000 

-  .  6  X  *  * 

•  w  V  - 

•  C  ;  3 

.ooor 

-..*:2 

.  220  * 

.12* 

-.1671 

-.2199 

•  6  e  «•  -* 

*  •  C  7 «  a 

- .  25  6  ? 

,  *  l  t  : 

.:  7f 

-.*097 

-..’92 

.1*1: 

.353 

-.61x1 

-.2211 

.  326  : 

.x;  = 

-.60*3 

-.3227 

.  2  29  * 

.x  36* 
.6?x’ 
.6129 
.6X0'. 
.2539 
.16  9' 

.  l’C? 
• .  1  i  x 
-.169’ 
-.111* 

.oro  : 


.oroo 
.  oco: 

2.7«C3 
2.721c 
2.7260 
2.9SX-. 
3.2X01 
3.X-6  5 
x.xP55 
5.X  113 
6.2X21 
6.5654 
5. 1X07 
X.  X  1  9  9 

3.056  1 
2.9.24 
3.  0  72  ! 

•  oco: 


L  •  x  •  2 . 

U*3 

10  OCO 


- .  L06  5 
.0*06 
-.0*27 

-.0327 

-.,529 

-:x?c£ 

-.1277 
-.2256 
- . 1 6  66 
-  •  L°  1  ’ 
-.0X37 
-.0276 
■.ul£ 7 

-  n  .  . 

-Iclh 

.0000 

:tfJ2 


3.3662 

2.756’ 

2.6161 

2 . ex  ol 

2.6966 

f I?ett 

3.09x3 

3. X665 
5.101" 
7.529* 

IXti 

.0000 

12.6572 

..002 

2.98x6 

.0002 


Table  B12B. 


l  • :  x  *  / 
u  r  7 

1CC0-  > 


.3’:- 
.5  7s. 
•  6^.* 
i  .x:.. 

1  .994. 

2  .7 7s 7 
2  .66  1. 
2  .6917 
1  .7905 

.2  C  3  x 


•c::o 

•  u  .  -J  w 

3  .1  6t  1 
3  ,30s- 

3 .1 9xx 

3  .2  x  .1 
3  ,x  *  2o 

x .o*!: 

5  .5  6  s  7 
X .661. 
5  .317. 
x .6730 
X  ,X  7eo 
3  .92 .s 
3  .l*9t 
3.10  7c 
3 .0557 
.00,. 


u*  v • :/ 
or: 
io:  oo 


.00.. 

-.12:: 

-.ir:-. 

-.lot: 

:Xn 

-.1  ’:x 
-.1677 
-  .  0  ’  t  5 
-.O’ct 
-.01x7 
-,o':7 

-.0*6  5 


»  »  -  *  - 
.  *  X  9  . 
>tx.« 
.  6 1 : 4 

.X  X  :  5 

,321c 
.  2  ■* :  2 


3.o*-;x 
3 . 0  7s  x 

3.1*15 
3.3xt e 

3 . 2  *x  1 

3.3  73  . 
3 .5 "9  J 
3.7*56 
3.9X7’ 
3.9  196 
3.76x0 

.020  ; 
9.3-91 
.0*0  0 
2.99X6 
.0020 


Mean  Profile  Data 


x  =  68  in. ,  Te  =  1.6% 


r 

r/ 

U 

T 

L-UE 

N 

XNCHt  1 

Oil  1  £ 

87/ SEC 

roc.f 

L/Of 

TMtTA 

L'T  40 

U  l  *  ) 

V  (  *  I 

1 

•  CO  7  3 

.007 

41.32 

66.  26 

.421 

.371 

-14 .{ 98 

10.616 

7.672 

13.250 

2 

.0083 

.007 

43. 36 

85.45 

.442 

.403 

-14.352 

11.364 

6.544 

15.057 

3 

.0054 

.308 

45.51 

64. 6J 

.464 

.436 

-13.794 

11.922 

12.403 

9.256 

it.  r*.fc 

<4 

.0106 

.00  5 

47.  35 

84.65 

.4  82 

.426 

-13.313 

9.047 

19.215 

5 

.0113 

.310 

46.06 

64.49 

.490 

.442 

-13.123 

12.593 

9.387 

20.460 

6 

.0132 

.012 

50  •  CO 

8  T.  36 

.535 

.485 

-12.617 

13.099 

10.261 

23.«14 

7 

.0146 

.013 

51. 06 

8  3.24 

.520 

.485 

-12.334 

13.382 

10.360 

26.445 

8 

.015  3 

.014 

51.66 

83.27 

.526 

.486 

-12.179 

13.537 

10.35b 

2  7  .  7  1  C 

9 

.0176 

.016 

52.95 

82.67 

.535 

.504 

-11.645 

13.871 

10.685 

3  1  .  8 1  6 

10 

.0194 

.017 

53.42 

82.37 

.544 

.523 

-11.722 

13.951. 

11.103 

35.121 

11 

.0214 

.019 

54.41 

6  2.09 

.554 

.536 

-1 1 .462 

14.254 

11.380 

36.71b 

12 

.0232 

.021 

55.16 

81.52 

.562 

.541 

-11.266 

14.450 

11.477 

41.540 

13 

.0245 

.322 

55.70 

81.51 

.567 

.541 

-11.124 

14.552 

1 1 .4eS 

45.063 

14 

.0266 

.024 

56.  15 

81.83 

.572 

.  544 

-1 1  .007 

14.709 

11.553 

46  .  u97 

15 

.0285 

.325 

56.82 

81.71 

.575 

.549 

-10.632 

14.864 

11.652 

11.793 

51.570 

lb 

.0304 

.027 

57.06 

81.54 

.  5  B  1 

.556 

-10.769 

14.547 

55.005 

17 

.0322 

.025 

57.40 

81.34 

.565 

.564 

-10.680 

15.036 

11.957 

56.258 

IB 

.0366 

.035 

59. 25 

80.24 

.604 

.607 

-10.164 

15.532 

12.874 

70.18? 

19 

.0454 

.041 

60. 37 

60.03 

.615 

.615 

-9.901 

15.815 

1  3 . 04  S 

€2.11? 

20 

.0524 

.047 

61.65 

79.94 

.626 

.618 

-9 .568 

16.14o 

13.117 

94.772 

21 

•  OS  6  6 

.052 

62. 44 

79.73 

.636 

.626 

-9.361 

16.355 

13.293 

105.979 

22 

.0657 

.055 

63.60 

79. C2 

.646 

.654 

-9.056 

16.660 

13.880 

1 1  6  .  ?  1  2 

23 

.0725 

.065 

64.57 

79.39 

.656 

.640 

-6.803 

16.513 

13. sec 

131.104 

2<» 

.0787 

.070 

65.  19 

78.99 

.664 

.656 

-8  .640 

17.076 

13.911 

142.311 

25 

.0854 

.076 

66. 02 

78.53 

.673 

.673 

-8.422 

17.294 

14.267 

154.422 

26 

.0925 

.083 

66.63 

7e.  18 

.676 

.687 

-6.270 

17.44b 

14.579 

167.256 

27 

.0986 

.086 

67.23 

78.53 

.685 

.673 

-8.104 

17.612 

14.286 

1  76  .  282 

28 

.1053 

.094 

67.83 

77.97 

.691 

.695 

-7.946 

1 7 . 7b  8 

14.756 

190.393 

29 

.1124 

.103 

68.49 

77.70 

.6  56 

.706 

-7  .77b 

17.940 

14.979 

203.027 

30 

.116  6 

.106 

68.96 

77.44 

.703 

.716 

-7.645 

16.071 

15.192 

214. 4»4 

31 

.1255 

.112 

69.24 

77,70 

.705 

.706 

-7.579 

18.137 

18.361 

14.983 

226.906 

32 

.1327 

.118 

70.17 

77.71 

.715 

.706 

-7.335 

14.971 

239.920 

33 

.1494 

.133 

71.33 

76.75 

.727 

.743 

-7.032 

16.684 

15.769 

270.107 

34 

.1672 

.149 

72.67 

76.  74 

.743 

.743 

-6.679 

19.037 

15.773 

302.26? 

35 

.184$ 

.165 

73.35 

76.45 

.747 

.755 

-6  .501 

15.215 

16.016 

333.56! 

36 

.2026 

.181 

74.45 

76.11 

.756 

.768 

-6  .2  j  3 
-5.886 

19.503 

15.630 

16.302 

366.270 

37 

.2154 

.196 

75.73 

76.25 

.771 

.762 

16.178 

396.63? 

38 

.2377 

.212 

76.13 

7C.  89 

.775 

.777 

-5.773 

19.943 

16.479 

429.717 

39 

.2544 

.227 

77.  lb 

75.43 

.766 

.795 

-5  .498 

20.216 

16.663 

459.40! 

40 

.2725 

.243 

78.23 

75.43 

.757 

.795 

-5.232 

20.464 

16.659 

492.621 

41 

.2854 

lift 

78.68 

75.67 

.601 

.  785 

-5 .106 

20.610 

16.667 

523.169 

42 

.3076 

79.8b 

75.59 

.614 

.789 

-4 .795 

20.921 

16.732 

55b. 067 

43 

.3675 

.328 

8  2.06 

74.49 

.636 

.631 

-4  .215 

21.501 

17.642 

664 . 34  1 

44 

.4275 

•  3e2 

64.49 

73.71 

.661 

.862 

-3.584 

22.13 2 

18.291 

772.74b 

45 

.4875 

.435 

66 . 40 

73.22 

.883 

.881 

-3  .084 

22.632 

18.696 

881.251 

46 

.5476 

.4e5 

68.40 

77.09 

.901 

.886 

-2.558 

23.156 

18.804 

989.667 

47 

.6074 

.542 

90.  14 

72.73 

.918 

.901 

-2.103 

23.613 

19.125 

1097. 9£0i 

46 

.6675 

.596 

91.82 

72.70 

.535 

.901 

-1  .664 

24.052 

19.126 

1206.6 lb 

49 

.7274 

.649 

93. 24 

72.28 

.950 

.918 

-1  .292 

24.424 

19.476 

1314.690 

50 

.7876 

.703 

94.42 

71.77 

•  9  b2 

.938 

-.962 

24.734 

19.901 

1423.7Cfc 

51 

.84  74 

.756 

95. 34 

71.44 

.971 

.951 

-.741 

24.975 

20.176 

1531.CC0 

52 

.907  ; 

.810 

VS.  93 

70.79 

.977 

.976 

-.587 

25.129 

20.7C9 

1640.074 

S3 

.9676 

•  864 

96. 71 

71.17 

.985 

.961 

-.378 

11:111 

20.395 

1749.071 

54 

1.0276 

.917 

97.21 

7C.78 

.990 

.976 

-.251 

20.719 

1857.626 

55 

1.0876 

.971 

97.54 

7r.  70 

.994 

.980 

-.164 

25.552 

20.792 

1965.96  1 

56 

1.1476 

1.024 

97.95 

71.01 

.996 

.967 

-.055 

25.658 

20.529 

2074. -3b 

57 

1.2073 

1.078 

97.92 

70.27 

.997 

.996 

-.065 

25.651 

21.142 

2182. 349 

58 

1.2676 

1.131 

98.05 

70.31 

.995 

.995 

-.032 

25.664 

21.109 

2291  . ?  4  b 

59 

1.3272 

1.185 

98.13 

70.48 

l.CCC 

.986 

-.011 

25.705 

20.967 

2399.078 

60 

1.3874 

1.236 

98.20 

70.16 

1.000 

1.000 

.008 

25.724 

21.216 

2507.694 

61 

1 .4474 

1.252 

98.16 

70.54 

1.000 

.988 

.004 

25.720 

20.921 

2616. 349 

62 

1.SC7S 

1.345 

98. 25 

70.74 

1.0C1 

.978 

.024 

.033 

25.736 

20.753 

2 724.985 

63 

I'Mtt 

1.9C3 

98.30 

70.51 

1.001 

.987 

25.749 

20.942 

3854.724 

64 

2.461 

98.20 

70.  36 

1.000 

.992 

.007 

25.723 

21.055 

4984.644 

65 

3.3825 

3.015 

96.  16 

70.14 

1.000 

1.002 

.003 

25.719 

21.252 

6114.203 

66 

4.00  7  4 

3.577 

98.27 

70.59 

1.0C1 

.984 

.026 

25.742 

20.87e 

7243.761 

Table  B13 


RUN  NO.  3.  POINT  2. 

BOUNDARY  LAYER  PROPERTIES  STANDARD 

LINEAR  SUBLAYER 

INTERPOLATION  FUNCTION  FRC 
TO  VALL  -ALL  TO  Y*  =  2 

FREE  STREAM  VELOCITY  r  98.170  96.170 

FREE  STREAM  TEMPERATURE  =  70.18C 

HALL  TEMPERATURE  =  95.750 

•(ALL  HEAT  FLUX  =  .  07954 

FREE  STREAM  DENSITY  =  .07552 

FREE  STREAM  KINEMATIC  VISCOSITY  =  .T001619 

DENSITY  OF  FLUID  AT  WALL  =  .07204 

KINEMATIC  VISCOSITY  OF  FLUID  AT  WALL  =  .0001750 

WALL/FREE  STREAM  DENSITY  RATIO  =  .95396 

LOCATION  REYNOLDS  NUMBER  (REX)  =  3435559.84 

INPUT  VALUE  OF  VELOCITY  DELTA  :  1 . 3X30  OC 

INPUT  VALUE  OF  TEMPERATURE  DELTA  =  1.30000 

CALCULATED  DELTA  =  1.12046 

DELTA  99. 5X  INPUT  =  .0C30C 

DISPLACEMENT  THICKNESS  (DELSTAR)  =  .15429  .15443 

momentum  thickness  (theta)  =  ,1C938  .iQ9t7 

£Nt°GY-DI SSTPAT  ION  THICKNESS  =  .  19539  .  19556 

ENTHALPY  THICKNESS  =  .00597  .00598 

SHAPE  FACTOR  1?  (OELSTAR/THETA  )  =  1.41057  1.40619 

SHAPE  FACTOR  32  l  ENERGY /THET A  )  =  1.78634  1.76327 

momentum,  thickness  Reynolds  number  =  5526.16  5540.06 

DISPLACEMENT  THICKNESS  PEYNCLDS  NUMBER  =  7795.01  7o02.27 

SKIN  FRICTION  COEFFICIENT  =  .002865 

FRICTION  VELOCITY  =  3.81748 

LA.  OF  THE  WALL  CONSTANT  (K)  r  .41000 

L  A «  OF  The  wall  CONSTANT  (C)  =  5.CA0CU0 

-AKE  STRENGTH  =  .44004 

CLAUSERS  'CELTA*  INTEGPAL  =  -3.62559  -3.91767 

CLAUSEPS  »G»  INTEGRAL  =  24.97566  24.86643 

DISPLACEMENT  THICKNESS  -  CONSTANT  DENSITY  r  .14465  .14645 

MOMENTUM  THICKNESS  -  CONSTANT  DENSITY  r  .11055  .11085 

SHAPE  FACTOR  12  -  CONSTANT  DENSITY  =  1.30848  1.33924 


LOCATION  -X-  68.00000 

Te  =  1.6% 


Table  B14 


Fluctuating  Profile  Data 
x  =  68  in.,  Te  =  1.6% 


Y : 

N  INCHES 


4  •  DOC. 

:Si: 

.157 

.26; 

.j?: 

.400 
.520 
.600 
.75 : 


11 

.530 

1* 

1.050 

1’ 

1.2CC 

1* 

1  .35  J 

if 

l.soS 

16 

1.750 

17 

4  •  CC  u 

1* 

*n* 

•  4.-4- 

Y/ 

DO.  7 A 

j.  sc: 

•  Csu 

•  105 

•  16  6 
•  J  £  7 
.269 
.35  0 
.4*7 
.525 
.656 
.767 
.519 

1.C50 
1  .If  1 
1 .312 
1.531 
3.570 


.0137 

•  07  56 
.  y  7  5  1 
.0736 
.0677 
,  ufe  7" 

•  u  6  1 7 
.0612 
,C25<* 
.0«3; 
.0312 
.0225 
.our 

.0197 


.0125 

.cc:r 


.0179 
..755 
.  036? 
..379 
.0371 
.:75i 

.0399 

.0319 

.0000 

..?95 

.0152 

.0160 

.0163 

.0190 

.0191 

.0190 

.0163 


.0153 

.059  1 
.0962 
.0967 
.099  7 

.0937 

:8??; 

.0390 

•  C?00 
.029- 
.3156 
.  3165 
.015  r 
.0195 

•  0!9u 
.019  2 
.0207 


VJ  *  V « /UE 


•  0050 
.0363 
.0371 
.0335 
.0323 
.0329 
.0311 
.0765 
.0002 
.0202 
.0135 
.0066 
.0279 
.0029 
.0036 
.0096 
.0096 
.0002 


0  ’V  •  / 

6 ?»4. 


L*  V*  /!■•» 

-.0®c0 
•  3  5  7  «* 
.9  76  -7 
.3657 
.  3  7t  2 
.“j?- 
.4*33 
.36  7  0 
.00 
.399  7 
.3176 
.i'i: 

.1976 
.075  1 
*,2677 
-.1267 
- .0°59 


n  inc«e: 


Y/ 

DOLl  < 


v  *  t  *  /  C  7 


6 '7/6  2 


r:/t? 


9*2/72 


•  U  .  W 


3.627 

.299 

.125 


5 

.60  0 

.5U 

i" 

.75; 

.6  =  6 

*i 

.52: 

.  7  c  7 

i? 

1.050 

.5  19 

1? 

1.200 

1 .05  7 

i“ 

1.350 

l.lel 

ir 

1.500 

1.312 

i£ 

1.750 

1  .'31 

i* 

4  «  00  - 

i.i-z 

!• 

.000 

.cco 

-.0336 
.1162 
.1955 
*12  62 
•  U-6 
.  I3c« 
.1923 
.1151 
.  002C 
.  i?ir 

.095' 

.0637 

.0635 

.OCf 

-.0229 

-.C«25 

~:Soot 


.2537 
.6216 
.5567 
«  cC-9  2 
.f7cj 

.5c?t 

.5759 

.  6995 

.35  £3 
.5599 
•  5ri9 
.9909 
.  3f  91 
.3259 
.25  79 
.2753 
.25  13 
.2202 


.9316 

:\\l\ 
.1662 
.1722 
.162  = 
.1*46 
.1665 
.  0203 
.1777 
.1509 
.2297 
.2772 
.  2« 69 
.  3909 
•  3  5C  c 
.  925  1 
.0302 


.3196 

.2786 

.296? 

.2355 

.2525 

.2995 

.2962 

.29e7 

.6917 

.2679 

.3063 

.3355 

.3367 

.3752 

.3622 

.3659 

.3196 

.0000 


.0  395 

~  c  : 

•  1  -  * 

.0“1  3 

.m: 

.  u  5  : 
.2*13 
.2?:  ; 


::Pi- 


INCHES  0£lU  75T*/V*^*  7  *  /  TT  AU  TVITW-TEj 


1  9,300  3.500 

t  .252  .;?9 

9  .12-  .1-6 

9  .15  7,  ,U6 
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Mean  Profile  Data 
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RUN  no.  3.  POINT 

BOUNDARY  LAYER  PROPERTIES 


FREE  STREAM  VELOCITY 
FREE  STREAK  TEMPERATURE 

wall  temperature 

•ALL  HEAT  FLUX 
FREE  STREAM  DENSITY 
FREE  STREAM  KINEMATIC  VISCOSITY 
DENSITY  OF  FLUID  AT  WALL 
KINEMATIC  VISCOSITY  OF  FLUID  AT  WALL 
WALL/FREE  stream  density  ratio 
L  CCA  T ION  REYNOLOS  NUMBER  (REX) 
INPUT  VALUE  OF  VELOCITY  DELTA 
INPUT  VALUE  OF  TEMPERATURE  DELTA 
CALCULATED  DELTA 
DELTA  99. 5X  INPUT 

displacement  thickness  (oelstar) 
momentu"  thickness  (Theta) 

£N£PGY~ DISSIPATION  THICKNESS 

enthalpy  thickness 
shape  factor  12  (delst ar/thet a) 
shape  FACTOR  32  ( ENERGY /THET A  ) 
MOMENTU"  THICKNESS  REYNOLDS  NUMBER 

displacement  thickness  Reynolds  number 

SKIN  FRICTION  COEFFICIENT 
FRICTION  VELOCITY 
LA.  OF  THE  WALL  CONSTANT  <K> 
LAw  GF  THE  WALL  CONSTANT  (C) 
WAKE  STRENGTH 

CLAUSERS  'DELTA’  INTEGRAL 
CLAUSERS  'G'  INTEGRAL 
DISPLACEMENT  THICKNESS  -  CONSTANT  DENSITY 
MOMENTUM  THICKNESS  -  CONSTANT  DENSITY 
SHAPE  FACTOR  12  -  CONSTANT  DENSITY 


LOCATION  -X- 
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Figure  B12E.  Boundary  Layer  Flatness  Distributions 
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Figure  B14E.  Boundary  Layer  Flatness  Distributions  x  =  52  in,  Te  =  4.7 Z 
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Figure  B15.  Mean  Velocity  and  Temperature  Profiles  x  *  68  in.,  Te  *  4.2* 
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UV  TRIPLE  PRODUCTS 


Figure  816C.  Boundary  Layer  Triple  Product  Distributions 


TRIPLE  CORRELATION 


Figure  B16D.  Boundary  Layer  Skewness  and  Triple  Product  Correlation 
Coefficient  Distributions  x  =  68  in,  T„  = 
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Table  B25 


RUN  NO.  1.  POINT  2. 

BOLNDARY  LAYER  PROPERTIES  STANOARD 

LINEAP  SUBLAYER 

interpolation  FUNCTION  from 
TO  W  LL  w  ALL  TO  Y*  =  35 

FREF  STREAM  VELOCITY  =  98*967  96.967 

FREE  STREAM  TEMPERATURE  =  71.450 

HALL  TEMPERATURE  =  95.270 

HALL  HEAT  FLUX  =  .07956 

FREE  STREAM  DENSITY  =  .07504 

FREE  STREAM  KINEmATIC  VISCOSITY  =  .0001633 

DENSITY  OF  FLUID  AT  wALL  =  .07182 

KINEMATIC  VISCOSITY  OF  FLUID  AT  HALL  =  .C001764 

*  ALL/FREE  STREAM  CENSITY  RATIO  =  .957J8 

LOCATION  REYNOLDS  NUMBER  (REX)  =  424417,8.31 

INPUT  VALUE  OF  VELOCITY  DELTA  =  2.05000 

INPLT  VALUE  OF  TEMPERATURE  DELTA  =  2.10D0C 

CALCULATED  OELTA  =  1.56079 

DELTA  99.5*  INPUT  =  .00000 

displacement  thickness  (oelstar)  =  .17986  .16012 

MOMENTUM  THICKNESS  (THETA)  =  .13446  .13465 

ENERGY-OISSIPATION  THICKNESS  =  .24493  .24496 

enthalpy  thickness  r  .00715  .C07is 

SHAPE  FACTOR  1?  ( OEL ST AR/Th£T A  )  =  1.33749  1.33772 

SFAPE  FACTOR  32  ( ENERGY /THETA  )  r  1.82137  1.61943 

momentum  thickness  reynclos  number  =  6794.61  6803.22 

DISPLACEMENT  THICKNESS  REYNOLDS  NUMBER  =  9087.75  9100.76 

SKIN  FRICTION  COEFFICIENT  =  .002989 

FFICTION  VELOCITY  r  3.91146 

LAW  OF  THE  WALL  CONSTANT  (K)  =  .41D00 

LA  H  OF  THE  WALL  CONSTANT  (C>  =  5.00000 

HAKE  STRENGTH  r  .17313 

CLAUSERS  ’DELTA’  INTEGRAL  =  -4.22327  -4.37735 

CLAUSERS  ’G’  INTEGRAL  =  23.76286  23.61295 

DISPLACEMENT  THICKNESS  -  CONSTANT  DENSITY  =  .1698C  .17297 

MOMENTUM  THICKNESS  -  CONSTANT  DENSITY  =  .13561  .13579 

SHAPE  FACTOR  12  -  CONSTANT  DENSITY  =  1.2521?  1.27363 


LOCATION  -X-  64. OOOOC 

Te  =  3.9% 


Table  B26 


Fluctuating  Profile  Data 
x  =  84  in. ,  Te  =  3.9% 
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Fluctuating  Profile  Data 
x  =  84  in.,  Te  =  3.9% 
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